BAND-EDGE SOLITONS, NONLINEAR SCHRODINGER / GROSS-PITAEVSKII 
EQUATIONS AND EFFECTIVE MEDIA 

B. ILAN* AND M. I. WEINSTEINt 

Abstract. We consider a class of nonlinear Schrodinger / Gross-Pitaevskii (NLS/GP) equations with periodic 
potentials, having an even symmetry. We construct "solitons", centered about any point of symmetry of the poten- 
tial. For focusing (attractive) nonlinearities, these solutions bifurcate from the zero state at the lowest band edge 
frequency, into the semi-infinite spectral gap. Our results extend to bifurcations into finite spectral gaps, for focusing 
or defocusing (repulsive) nonlinearities under more restrictive hypotheses. 

Soliton nonlinear bound states with frequencies near a band edge are well-approximated by a slowly decay- 
ing solution of a homogenized NLS/GP equation, with constant homogenized effective mass tensor and effective 
nonlinear coupling coefficient, modulated by a Bloch state. 

For the critical NLS equation with a periodic potential, e.g. the cubic two dimensional NLS/GP with a periodic 
potential, our results imply: 

• The limiting soliton squared L 2 norm, as the spectral band edge frequency is approached, is equal to V e dge = 
C* x V C r, where V C r denotes the minimal mass soliton of the translation invariant critical NLS. V C r is also known 
as the Townes critical power for self-focusing of optical beams. 

• The constant is expressible in terms of the band edge Bloch eigenfunction and the determinant of the effective 
mass tensor, f the potential is non-constant, then < < 1 and V^ge. is strictly less than, V C r- 

The results are confirmed by numerical computation of bound states with frequencies near the spectral band edge. 
Finally, these results have implications for the control of nonlinear waves using periodic structures. 
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1. Introduction and Outline. "Solitons" are spatially localized concentrations of en- 
ergy, which are of great interest in many nonlinear wave systems. They arise from a balance 
of dispersion (or diffraction), which tends to spread energy and (focusing / attractive) non- 
linearity which tends to concentrate energy. Although their importance was first recognized 
in the context of hydrodynamics [66, 2], soli ton-like coherent structures are now understood 
to play a central role in contexts ranging from optical pulses (temporal solitons) to stationary 
beams (spatial solitons) of nonlinear optics [11, 40] to soliton matter waves in macroscopic 
quantum systems [46]. Advances in the design of micro- or nano-structured media have 
greatly enabled the control of optical and matter waves. Thus it is of interest to develop a 
fundamental understanding of the effect of inhomogeneities in a medium on the dynamics of 
nonlinear dispersive waves and, in particular, on the dynamics of solitons. See, for example, 
[24] for an experimental investigation of solitons in periodic structures. 

In this article we consider solitons in non-homogeneous media governed by a class of 
nonlinear Schrodinger / Gross-Pitaevskii (NLS / GP) equations: 

id t tp = -Atp + V(x)V> - H 2<T V> . (1.1) 

Here tp = ip( x >t) denotes a complex-valued function of (x, t) € R d x R 1 , d > 1. The 
potential V(x) is assumed to be a real-valued, smooth, periodic and symmetric about one or 
more points 1 2 
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1 Periodic potentials are often called "lattice" potentials. 

2 The main results of this paper extend to general nonlinearities of the form 

f(M 2 )ip = g[l + 0(W 2 )] M 2 °ip. 
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NLS / GP is a Hamiltonian system, expressible as: 



■a i 5H 

idttjj 



-HU/,,1/,*] = /" w • W + V(x)v^* l — (V>T + \ (1.2) 

7 (7 + 1 

where -0* denotes the complex conjugate of tp. By Noether's Theorem, the invariance t i-> 
< + to implies the time-invariance of H for solutions of NLS / GP. Furthermore, the invariance 
ijj h-> e' e ?/> implies the additional time-invariant quantity: 

P^V*] = f i>^* rfx. (1.3) 



The parameter a > allows for variation of the strength of the nonlinearity. In physical 
systems, we typically have a = 1. Allowing cr to vary enables one to quantify the balance 
between nonlinear effects and dispersive / diffractive effects, which depend on spatial dimen- 
sionality, d. Local well-posedness in time for the initial value problem for (1.1) with data 

V(x,t = 0) = o (x) (1.4) 

in Vo G H 1 ^) (see, for example, [16, 59]) holds for all cr > for d = 1,2 and all 
< cr < 2(d — 2)" 1 for c? > 3. Global well-posedness for arbitrary data holds for a < 2/d. 
For well-posedness for data in spaces of weaker regularity, see [10, 60]. 

The NLS/GP equation, (1.1), with a = 1 governs the dynamics of a macroscopic quan- 
tum state, a Bose-Einstein condensate, comprised of a large collection of interacting bosons 
in the mean-field limit [46, 25, 15]. The attractive nonlinear potential, — 1?/'| 2 , corresponds to 
a species of bosons, whose two-particle interactions have a negative scattering length. A sec- 
ond important area of application of NLS/GP is its description of the evolution of the slowly 
varying envelope of a stationary and nearly monochromatic laser beam propagating through 
a nonlinear medium [40, 11]. Here, the attractive nonlinear potential is due to the Kerr non- 
linear effect; regions of higher electric field intensity have a higher refractive index. In this 
setting, t denotes the distance along the direction of propagation and x G R 2 the transverse 
dimensions. In the quantum mechanical setting the potential, V(x), is determined by mag- 
netic and optical effects that are used to confine a cloud of bosons. In optics, the potential is 
determined by the spatial variations of the background linear refractive index of the medium. 
The functional V denotes the optical power or, in the quantum setting, the particle number. 

Nonlinear bound states or solitons of NLS/GP are solutions of standing wave type: 

^(x^-e-^x^) (1.5) 

where /i denotes the frequency (propagation constant in optics, chemical potential in quantum 
many-body theory) and u is a real-valued solution of 



(-A x + V(x) ) u(x, fx) - u 2CT+1 (x, M ) = M u(x,/x), u(;h) e H^W 1 ). 
We shall construct solutions of (1.6) with fj, located in a spectral gap of — A + V. 



(1.6) 



The properties of solitons in homogeneous media, V = 0, are reviewed in detail in Sec- 
tion §2. Briefly, for cr < 2/d (subcritical nonlinearities) dynamically stable solitons exist at 
any prescribed L 2 norm (in one to one correspondence with any fj, < 0). In the critical case 
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and supercritical cases, a > 2/d, soli tons are unstable. 

We raise several motivating questions and outline our results in the next subsection. 

(Ql) Persistence and stability: What is the effect of a periodic potential on the existence 
and stability properties of solitons? 

(Q2) Stabilization: Can unstable solitons be stabilized by a potential, V^(x)? This ques- 
tion was first addressed in [48] in the context of localized potentials and more re- 
cently for more general potentials [54], e.g. periodic, quasiperiodic. 

(Q3) Excitation thresholds / Minimal mass solitons: How does periodic structure effect 
soliton excitation thresholds? For critical nonlinearity, a = 2/d, and V = the 
soliton squared L 2 norm is independent of //; 'P[it(-,/i)] = V C r- Thus, there is 
an L 2 threshold below which there are no solitons. This L? excitation threshold 
for soliton formation is of great physical interest [65, 26, 19, 41, 9]. In optics it 
corresponds to the critical power for self-focusing [18, 42]. Such solitons are also 
often called minimal mass solitons. See Remark 3.3 . 

1.1. Outline of Results. In order to outline the results of this paper, we begin with very 
quick review of the spectral theory of Schrodinger operators, — A + V, for V periodic [47, 
23, 36]. If V(x.) is a periodic potential, then the spectrum of — A + V is real, bounded below, 
tends to positive infinity, is absolutely continuous and consists of the union of closed intervals 
( spectral bands). The open intervals separating the spectral bands are called spectral gaps. 
One dimensional Schrodinger operators with periodic potentials generically have infinitely 
many gaps. In dimensions d > 2, there are at most finitely many gaps. 

We denote by E* lowest point in the spectrum, the left endpoint or edge of the first 
spectral band. E* is simple and is the ground state (lowest) eigenvalue of — A + V, subject to 
periodic boundary conditions on the basic period cell of V. The eigenspace associated with 
E* is spanned by w(x), a non-trivial solution of: 

(—A + V(x)) u>(x) = E*w(x), w(x) periodic 

For the case V = 0, -E 1 * = and we can take u>(x) = 1. 

The present work considers the bifurcation and dynamic stability properties of families 
of solitons emerging from a spectral band edge. Such edge-bifurcating solitons have a multi- 
scale character described below. Our results include the following: 

1. Theorem 3.1: Let xo denote any point of symmetry of V(x) 3 . 
There is a family 

/j i y u(x, (i ) w (E* — p)^ F ( y/E* - fj, (x - x ) ) u>(x), 

which bifurcates from the zero solution at energy E* into the semi-infinite gap 
(— oOjE*) for < E+ — fi sufficiently small. Here, a G N for d = 1,2 and 
(7 = 1 for d = 3. 

F(y) denotes the soliton profile for an effective medium with effective mass tensor, 
A % i given by (3.5) and effective nonlinear coupling constant, 7 e ff given in (3.6), and 



3 i.e. /(x) = f(xi, ... ,x d ) is symmetric (about the origin) if f(x lt ... ,x d ) = f(s 1 x 1 ,...,s d x d ), Sj = 
ztl. xo is a point of symmetry of V(x) if V(z) = V(xo + z) is symmetric. Thus, by translating coordinates, we 
can arrange for a point of symmetry to be at the origin. 
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satisfies the homogenized soliton equation: 



- d y ^d yj F(y) - 7eff F 2 °+\y) = -F(y) . (1.7) 

F > 0, Fei/ 1 (M d ) 

The leading order expansion is constructed via multiple scale expansion. The error 
term is studied by decomposition of the corrector into spectral components near and 
far from the band edge, and estimated via a Lyapunov-Schmidt strategy; see also [14, 
21, 22]. The results can be extended to solitons near edges of finite spectral bands 
for focusing and defocusing nonlinear potentials under more restrictive hypotheses 
on V; see Section §3.1. A variant of Theorem 3.1 holds in dimension one, in any 
spectral gap, near a "positive curvature" band edge; see Theorem 3.6. 
Corollary 3.3, part 1 : Consider the critical cases: a = 1 , d = 2 and a = 2, d = 1. 
Near the band edge, i.e. for E* — \i > and small, we have: 

V[u(-,n)] = C*V cr + (m-£*)Ci* + 0((m-£*) 2 )- (1.8) 

Here, 

(a) V C r = F[R(-, — 1)], where R(-, — 1) denotes the unique (up to translations) 
solution of 

AR - R + Ri +1 =0, R > 0, R e H 1 ; (1.9) 

see (Q3) above. 

(b) 

■Pedpe = lim P[lt(-,/Lt)] = C* Vcr 
fj,— >E* 

is given by (3.15) and satisfies the inequality < £* < 1, unless V is identi- 
cally constant. In the latter case, = 1. 

(c) 



Ci. = - 



V[u(;fl)] (1.10) 



is given by (3.16). 

For periodic potentials, V, of the form 6 V(x), where S is sufficiently small, 

we show that £i* > 0; see (3.18). 

Positive slope conjecture: In general, (i* > 0. 
Both and (i* depend on the edge (periodic) Bloch eigenstate and the Hessian 
matrix (of 2 nd partial derivatives) of the band dispersion function, D 2 Ei, near E*. 
The latter is often called the inverse effective mass tensor. 
3. Instability for \x near the band edge: Consider the critical cases a = 1, d = 2 and 
a = 2, d = 1. For V nonzero and /i near E*, nonlinear bound states are linearly 
exponentially unstable, provided £i* > 0. We conjecture d* > 0, in general, and 
have verified it for potentials V = SV, with S sufficiently small. 
For V = 0, the linear instability is algebraic, although for the nonlinear dynamics, 
solutions can blow up in finite time or decay to zero dispersively (diffractively) as t 
tends to infinity. In contrast, since for V ^ 0, for fj, close enough to the spectral band 
edge, the curve /j, H> V[u(-,fj,)] lies below the line V = Vcr, although a solution 
with data near a soliton, u(-,fi), with frequency near the band edge does not remain 
nearby in H 1 , the solution exists globally in time in H l . 
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4. Section §6: Numerical computations are used to illustrate the asymptotic results 
and to study the global behavior. 
Figure 1.1 below summarizes a key consequence of our results. 




-10 2 -10 1 -10° 

-log(-n) 

FIG. 1.1. Plot of power curves: V[u(-,fi)] vs. fi (using semi-log axis) for the quintic one-dimensional ' NLS/GP 
equation, (6.3), with Vo = 10 and K = 2tt (here E* Rj — 1 . 23 J. Solid (blue) curve corresponds to power curve 
for soliton family centered at a local minimum. Dashed (red) cun>e corresponds to centering at a local maximum. 
Agreement is shown between numerical computations and the analytically obtained value for the band edge power 
(dashed / black horizontal line), V e j_ ge = lim^^^ V[u(-,p)] = £ * X V cr ~ 2.2 [Eq. (3.14)]. For p, large and 
negative V(p) converges to Vcr ~ 2.72 (dashed /green horizontal line), which is the critical power of the Townes 
soliton in translation invariant (V = const) case. 

For any nontrivial periodic V(x) the limiting L 2 norm at the band edge is strictly less 
than that of the homogeneous medium. The slope of the curve, p — > V [u(-, p)} is strictly 
positive. As — p = \fi\ increases, solitons become increasingly localized in space, and thus 
depend more and more on the local properties of V. The limiting —> oo) squared L 2 
norm is V cr . The orbital stability theory, outlined in Section §2 implies that solitons with 
energies /j, near the band edge (where d^V \u{-,n)\ > 0) are unstable, while those which are 
centered and sufficiently concentrated (— ft sufficiently large, d^V [u(-, /i)] < 0) about a local 
minimum of V are stable. It is natural to conjecture that for localized initial conditions with 
L 2 norm strictly less than inf^<s, V [u(-, (J,)], solutions to the initial value problem disperse 
to zero as t — > oo; see the discussion in the proof of part 3 of Theorem 3.4 and [64]. 

REMARK 1.1. Concerning the dependence of y, i— > V[u{- , fj,)] for [i near the band edge, 
on the nonlinearity parameter, a, and dimensionality d (see Theorem 3.2 and Corollary 3.3) 
it is useful to recall the analogous behavior in the translation invariant case: V = 0. In this 
case, NLS is also invariant under dilation: 

ip(x,t) i — ^ \iip{\x,X 2 t) ■ (1.11) 

Let R{-, ft) denote the positive (unique up to translation), solution of 

-AR- R 2a+1 = fiR, 

By uniqueness 

i?(x,/i) = |/x|*i?(|/i|3 X ,-l) (1.12) 

It follows that 

V[R(;^] = \\R(;^)\\t = \^-i \\R(;-l)g 
implying that as p — > E t , V[R(-, p)\ tends 
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• to O,fora < 2/d 

• to \\R(-;—l)\\%for a = 2/d, and 

• to +00, for a > 2/d; 

see figure 1.2. In one space dimension, the family ofsolitons is given explicitly by: 

R(x,p.) = [0 + 1) |/i|pscch- (erv^ x \ . (1.13) 
In the critical case, a = 2, 

7W,A0] = ^T [ scc Hv) dy = tt ~ 2.7207 ; 

see figure 1.2. Theorem 3.2 implies a similar trichotomy of behaviors for states bifurcating 
from the band edge, E*, of a periodic potential. Also, for a = 2/d, the curves fi 1— > V[u(-, fx)] 
in figure 6.3, are seen to be deformations (for minimum and maximum centered solitons) of 
the horizontal line /1 1— > V C r for the case 1^ = 0. 




FIG. 1.2. (A) L°° norm and (B) squared L 2 norm (V) as functions of frequency, {1, for the ground state 
solution ofEq. (1.6) in one dimension (d = I) with three nonlinear exponents: subcritical (crd = 1 < 2), critical 
(ad = 2), and supercritical (crd = 3 > 2); see legend. In L°°, bifurcation appears from a state with zero norm at 
E„ = 0. In L 2 , the limiting behavior as p — > 0~ depends on crd. 



Previous work: Formal expansions and numerical approximation of nonlinear bound states 
near spectral band edges for periodic and aperiodic structures and their linearized stability 
properties were presented in [57, 7, 67, 45, 12, 8, 49, 50, 17, 1, 13, 5, 52, 56]. The band edge 
limit of V, for case of a 2-dimensional separable potential was obtained by formal perturba- 
tion theory and numerically in [51]. Two-scale convergence methods have been applied to 
rigorously derive homogenized effective equations, valid on large but finite time scales, in [6], 
for the linear Schrodinger equation, and in [55] for the time-dependent NLS/GP, with two- 
scale type initial conditions. Bifurcation of localized states from the continuous spectrum 
into spectral gaps has been considered in [37, 38, 30, 29, 4, 58, 43]. The connection with 
nonlinear coupled mode equations is explored in [14, 44, 21, 22]. The Lyapunov-Schmidt 
strategy applied herein is motivated by these latter approaches. 

Outline: The paper is structured as follows. In Section §2 we discuss background for the 
formulation of our results. We state our main results in Section §3. In Section §4 a formal 
homogenization / two-scale expansion of solitons with frequencies near the band edge is de- 
rived. The expansion and error estimates are proved in Section §5. In Section §4.3 we derive 
the consequences of our expansion of band-edge solitons for the character of 7 : '[w(-,/i)] as 
p, — >• E*. Section §6 contains a discussion of numerical simulations illustrating our main 
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theorems. Section §7 contains a short summary and discussion. The latter sections of the 
paper are appendices containing technical results on the effective mass tensor. 

1.2. Notation. 

1. We shall write 77(e) = C(e°°) if 77(e) = 0(e q ) for all q > 1. 

2. Fourier transform of G: G(k) = / e - 2 ™ k - x G(x) dx 

3- x(° < |k| < 6) = characteristic function of the set {k : a < |k| < b} 

4. x (|V y | <a)G = J e 2 ™ k y X (|k| < a) G(k) dk 

5. iP, Sobolev space of order s; H s even space of even iT s functions 

11/11*. = E W Qa ^ ~ H/lli».« 

|a|<s 

6. H s sym , symmetric 7J S functions, 

/.e. / e i? s s ym if / e H s and f(x l7 . . . ,x d ) = f(s 1 x l7 ...,s d x d ) 7 Sj = ±1. 

7- n/iii MB ) = ^ i/wi 2 (i + izi 2 rdz 

8. GJ™(]R d ), functions in G m (R d ) with limit equal to zero as |x| -> 00 

9. $ denotes the fundamental period cell and B* the dual fundamental cell, or Brillouin 
zone. 

2. Background. 

2.1. Solitons and Stability Theory. We give a very brief review of the stability theory 
of solitonsofNLS/GP, (1.1). 

DEFINITION 2.1. The nonlinear bound state u(x, /i) ofNLS/GP is orbitally stable if for 
all e > 0, there is a 5 > smc/i f/iaf if the initial condition tpo satisfies 

inf || Vo ~ u(-,M)e 47 ||#i < 5, 

7G[0,2t) 

f/ien f/ie corresponding solution, ip(-,t), satisfies 

inf HVK-.t) - u(-,^)e ll \\ H i <e, for all ^ 0. 

7G[0,2tt) 

This notion of soliton stability for NLS is natural since NLS/GP, for V non-constant is in- 
variant under the group of phase translations, i\> H> e l9 ip, but not spatial translations. 

A central role in the stability theory is played by the operator 

L+ = -A -fx+ V- (2cr + l)u 2a , (2.1) 

the real part of the linearization of NLS/GP about «(•, /j). Let n_ (£+) denote the number of 
negative eigenvalues of L + . If it(x, /i) is a nonlinear bound state with /i < E* (frequency ly- 
ing in the semi-infinite gap) then n_ (L + ) < 00 and the following nonlinear stability theorem 
holds [63,48,64,28,54] 
Theorem 2.2. 

1. Let u(x, fi) denote a positive soliton solution of NLS/GP with fi in the semi-infinite 
gap (— 00, E 1 *). The nonlinear bound state, -0(x, i) = u(x, /j)e _IM * is orbitally 
stable if the following two conditions hold: 
(a) Slope (VK) condition: 



— V[u{-,n)] <0, and 
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(b) Spectral condition: L + has no zero eigenvalues and 



n-(L+) = l. 



(2.2) 



2. If either V[u{-, fx)] > or n_ (£+) > 2 then the soliton is unstable (nonlinearly 
unstable as well as linearly exponentially unstable). 
REMARK 2.1. As discussed in [54], the spectral condition can be associated with the 
suppression of a drift instability and the slope condition with the suppression of an energy- 
concentrating self-focusing instability. 

2.2. Spectral theory for periodic potentials. We consider the Schrodinger operator 
— A + V(x) acting in L 2 (M. d ), where V(x.) is smooth, real-valued potential which is periodic. 
That is, V(x + q) = V(x) for all x € R d . Here q = {qi, . . . , q^} denotes a linearly 
independent set of vectors in M. d that spans (over the integers) a lattice denoted by T. The set 



3 = 1 



1 1 

2' 2 



(2.3) 



is called a fundamental period cell. The first Brillouin zone B* is generated by the dual basis 
(n, . . . , Y d ) given by Yj • q k = 2ir S jk , i.e. 



B" 




1 1 

2' 2 



and the dual lattice, V*, is the integer span of B*. 

It is useful to review some well-known results of Floquet-Bloch theory [35, 47, 23, 36]. 
The spectrum of — A + V, denoted a(— A + V), consists of a union of closed intervals called 
spectral bands separated by gaps (also known as band gaps and photonic band gaps). The 
spectral bands are characterized as follows. 

For each k £ B* we seek solutions of the linear eigenvalue problem 



of the form w(x; k) 



( -A + V(x) ) u = Eu 
e lk 'Xp(x; w h ere £,( x . k) is periodic in x: 

2 



(2.4) 



(V + iky + V(x) p(x; k) = E(k)p(x; k) , 

p(x + qj;k) = p(x;k), j = l,...,d. 

For each k E B* this periodic elliptic boundary value problem has a sequence of discrete 
eigenvalues or band dispersion functions tending to positive infinity: 

Ei (k) < E 2 (k) < ■ • • < E m (k) < As k varies over the Brillouin zone B* each E m (k) 

sweeps out a closed subinterval of the real axis. The spectrum of — A + V acting on L 2 (M. d ) 
is the union of these subintervals: 



U m >i {E m (k) :ke B*} C [rnin V,oo) 



a(-A + V) 

and the states {it m (x; k) = e Jx k p m (x; k)} are complete in the sense that 

feL 2 (R d ) /(x) = y;/ (u m (-;k)J(.)) L2m Um (x;k)dk (2.5) 



We denote the lowest point in the spectrum of Eq. (2.4) and corresponding periodic 
eigenstate by 

£*=£a(0), iu(x) = p(x;k = 0). 
E„ is simple. We will often make use of the relation 

L*w = , w > , w(x + qj) = w(qj) , (2.6) 

where 

L* = -A + V" - -E* . (2.7) 

Thus, w is the periodic ground state of L*, > and is a simple eigenvalue of L* with 
kernel spanned by u>. Note that if P- 1 is the orthogonal projection onto the subspace {to} , 



P ± 9 = 9 - (ttW ) 71 — H" ' (/>#) = / /( x )fl( x )^ x - 

\IHI / IMI Jb 



(2.8) 



then 1 P- 1 is bounded on the space of L 2 periodic functions with fundamental period cell 
B. 

Finally, note that we may, without loss of generality, restrict to the case where the fun- 
damental period cell is [— vr, Tr] d . Indeed, if B is the fundamental period cell (see (2.3)), then 
define the constant matrix Q to be the matrix whose j th column is (27r) _1 qj . Then, under 
the change of coordinates x n- z = Qx, we have 

— V x ■ V x + V(x) acting on L 2 er {B) transforms to 

- V z • a V z + V(z) = - a v + ^( z ) 

i,j=l 1 3 



acting on L 2 per ([— tt, 7r] d ) where 

-, V(z)^V(Q- 1 z) 1 x = 



QQ T 
\detQ\ 



T 1 * 



3. Main Results. In this section we state our main results on bifurcation of solitons 
from the band edge, E = E„, into the semi-infinite gap. 



Hypotheses: 

(HI) Potential: V(x) is smooth and periodic with B = [— ir, Tr] d . 
(H2) Dimension / Nonlinearity 4 : 

d= 1,2 : a G N, d = 3 : a = 1 

THEOREM 3.1. Let xo denote any point of symmetry ofV(x). 
1. For all \i less than and sufficiently near E*, there is a family of nonlinear bound 
states of NLS/GP ("solitons"), £ iJ s (R d ), s > d/2, which is centered at 

Xo. 

4 The assumption on the nonlinear term can be made less restrictive. However, since to some of our results 
concerning the higher order character of fj, i— f V[u(-, fi)] depends on the construction of a multiple scale expansion 
to a sufficiently high order, we require a certain degree of smoothness of the nonlinear term in a neighborhood of 
zero. Note also that the methods and our results extend easily to more general nonlinearities, e.g. IC[\u\ 2 ]u (local or 
nonlocal). 
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2. These solutions bifurcate from the zero solution at band edge frequency fj, = E„ into 
the semi-infinite gap. Specifically, this family is given by the two-scale expansion for 
small e, 

Me = ^-e 2 , (3.1) 

u e (x, /Lt e ) = e- [w(x)F (e(x - x )) + elli (x, e(x - x )) (3.2) 
+ e 2 [/ 2 (x, e(x-x )) + rj(x; e)] , 

where ?}(x; e) satisfies the estimate for any s > d/2 

N(-;e)lk. < C e 3 (33) 
77ie ferazs ;« f/ie expansion are given as follows: 

u>(x) is f/ie banc/ ec/ge Bloch state [see Eq. (2.6)] and F(y) is the ground state 
solution of the NLS equation in an effective medium: 

d 

- £ d yi A«d Vj F(y) - T eff F 2CT+1 (y) = -F(y) . (3.4) 

i J = l 

F>0, FeH 1 ^) 

The matrix A %] is the inverse effective mass tensor [35], expressible in terms of the 
band dispersion function, E\(\s.), as 

_ 4 (3 „,,L^) l|^ (k = 0); (3 . 5) 
(w,w) 2okiOkj 

see Appendix A. The effective nonlinear coupling constant is given by 

Teff = JB f 2( \ ' ■ (3.6) 
J B W ( x ) dx 

A 13 is a symmetric, positive definite constant matrix and its determinant, the product 
of inverse effective masses, denoted by 

— = det( A ij ) < 1, (3.7) 

with — 1 only ifV{~x) is identically constant; see [34] and Appendix B. 

3. F(y) is a rescaled ground state of the NLS equation (1.9) as 

F(y) = (—) fl(A-3Sy,-l) (3.8) 

VTeff/ 

where S is an orthogonal matrix that diagonalizes the effective mass tensor, i.e. 

S lk A kl Sij = Aij = diag(Ai, . . . , X d ) , (3.9) 

where Xi denote the eigenvalues of A 11 . 

4. Combining (3.1), (3.2) and (3.8) gives, for E* — p > and small: 

i 

u(x,p) = (J^fj [r(a-?S(Ap)? (x-x ),-l) U7(x) + 0(Afi)i 
Afi = E* - fi (3.10) 
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5. The 0(e) = O ^ Apj and 0(e 2 ) = O ( A/i) corrections are given, 
(using y = e (x — x ) and summation over repeated indices) by 

0(e) : L7i(x,y) = [d x ,w(^)} 9 Vi F(y) , (3.11) 

0(e 2 ) : C/ 2 (x;y) = C/ 2p (x, y) + w(x)F 2h (y) 

U 2p (x,y) = L' 1 [ ( % + ±d Xi L~ x d Xi ~ A«) to(x) ] %^F(y) 
+ L; 1 [^ 2CT+1 (x)- 7offW (x)] F 2ff+1 (y), 
if fiM.(y) = 5(y) , (3.12) 

where S(y) is given by 

S(y) = (w^w)- 1 [(w,(A y -l)U 2p (;y) + (2a + l)U^U 2p (;y)) 

+ a(2a + l) (wM^-'U^^y)) + 2{w, V x ■ V y [7 3 (-, y)>] (3.13) 

Theorem 3.1 is proved in Section §4.3. 

Using expansion (3.1) we can derive the asymptotic behavior for V{p) = V[u(-, p)} as p, 
E*. 

THEOREM 3.2. Let x denote a point of symmetry ofV and u(-, (l) a soliton given in 
Theorem 3.1. 

1. For p near the band edge V[u(-, p)] is given by 

V[u(;p)} 

= | n-Ek \ '~* 

x [C,V[R(- Cu(fi-E,) + 0((p-E*) 2 )] (3.14) 

where V[R(- , — 1)], the optical power of the homogeneous NLS ground state, de- 
pends on a and d: 



AR - R 2a+1 =-R, R>0, ReH 1 



and 



i V* ((/, 



III,, 



w 



° < c*^ ( — ) i T,„ 2 : +2 i < i. (3.i5) 



where -[ g = — —— [ g(x) dx, and the slope is given by 
J B vol(B) J 

Ci* = 4 E / [ d *M*)] I 2 dx ■ [ I d v, F (y) I 2 d y 

- £ w 2 (x) dx J ^F(y) + y • V y F(y)) S(y) dy, (3.16) 

where S(y) is given by Eq. (3.13). Note that to order (D(\E* — p\ x ) the expansion 
is independent of xq, the soliton centering. 

11 



2. Positive slope for small potentials: Let V(x) = 8Vi{x), where \8\ is sufficiently 
small and V\(x) is a smooth periodic function on K with a zero cell average. Then, 
in the critical case a = 2 



C*~l-8<5 2 / [(-d^Vi] 2 dx, (3.17) 

JB 

Cu ~ 34V3 n5 2 -f [(-d xxx )- 1 V 1 ] a dx . (3.18) 

JB 



Here, (— one/ (— 9 x;KC ) are respectively the second and third-order inte- 
gration operators in B acting on the space of zero average functions to itself. Hence, 
Ci* > for small potentials. 
3. Positive slope conjecture: £i* [V] > ifV is non-constant. 
Theorem 3.2 is proved in Section §4.3, except for part 2, concerning small potentials, 
which is proved in Appendix D. 

REMARK 3.1. Concerning equality in Eq. (3.15). When V(x) is constant then so is 

w(x). In that case E« = 0, £7 x (k) = k 2 , and m' 1 = dct |2 _1 D 2 .. k .Ei(0)\ = 1. 
Therefore = 1. 

REMARK 3.2. That^ < 1 can be seen by considering eachfactor in the definition (3.15) 
separately. First, by Holder's inequality the quotient in the second factor of (3.15) is bounded 
one with equality holding iff w = constant. Furthermore, w is identically constant if and 
only ifV = constant. Concerning the first factor in (3.15), by Theorem 3.1, < to" 1 < 1 
with equality holding ifV = constant. Therefore, < ^ < 1 with = 1 if and only if 
V = constant. 

In the critical case, an immediate consequence of Theorem 3.2 is the following result for 
critical nonlinearity (a = 2/d): 

COROLLARY 3.3. Consider the critical case a = 2/d; by hypotheses (H1)-(H2) this 
implies either (d, a) = (1, 2) or (d, a) = (2, 1). 
1. As \i — > we have 

V[u(-,n)] = C*V cr + Ci*(p-E*) + 0{iv-E*f) . (3.19) 

Here, V C r — "P\R{'i — !)]• Since ^ < 1 for any non-constant periodic potential, it 
follows that the limiting power at the band edge is strictly smaller than V cr , 

V edge EE lim V[U(-,H)} = C* Vcr < Vcr ■ (3.20) 

M — ^-E* 

Theorem 3.2 is proved in Section §4.3. The band-edge limiting behavior (3.20) is illustrated 

in Figure 1.1; see also Figure 6.3. 

Concerning the NLS / GP dynamics near solitons, we have the following: 

THEOREM 3.4. Consider the critical case a = 2/d; by hypotheses (H1)-(H2) this 

implies either (d,o~) = (1,2) or (d, cr) = (2,1). Then, if the positive slope conjecture of 

Theorem 3.2 holds, then 
1. 



dP[u(-,n)] 
d[i 



> (3.21) 



and it follows from Theorem 2.2 that for /i such that E* — fi > and sufficiently 
small, u(-, fi) is unstable. 
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2. In particular, for small periodic potentials, by Theorem 3.2, for /i such that E^—/i> 
and sufficiently small, u(-, fi) is unstable. 
To complement this information about stability / instability of solitons we remark on V cr 
and Vedge as they relate to well-posedness and blow-up / collapse. 

THEOREM 3.5. Denote by i?(x), the ground state ("Townes soliton" ) for V(yi) = 0. If 

V[M = y^o(x)| 2 dx< j i? 2 (x)dx = V cr (3.22) 

then solutions ofNLS/GP (1.1) exist globally in time; no singularity formation / no collapse. 

REMARK 3.3. Recall that in the spatially homogeneous case, V = 0, if in addition to 
(3.22) we impose the stronger assumptions: ipo G H 1 and |x|V>o G £ 2 , (fen ^(x, i) fe/icfc fo 
zero as t ^ oo for a range of p > 2 764/; see a/so [33] for scattering results in H 1 . 

V cr is thus called a soliton excitation threshold . Excitation thresholds also play a role in 
systems without critical scaling symmetry. See, for example, [65, 26] and [19, 41, 9]. 

For V non-zero, the picture which emerges from the above theorems and numerics (see, 
for example, figure 1.1) is quite different. The minimal mass (minimal power), band edge 
power and V = critical mass are related by: 

^min V edge. ^ V C r . 

Here, 

Vmi-a = V \u{- , fj. min )} = min :P[w(-,^)] , (3.23) 

fj,<E, 

where in (3.23): fj. i— > V[u(-, fi)] is computed along the family of solitons centered at a local 
minimum; see the solid curve in figure 1.1. Along this soliton curve, computations indicate 
thatu(-,fi) > and n_(L + ) = 1. By Theorem 2.2 (applied for V periodic) there is an open 
set of initial data in the phase space H : 

{V-0 G H 1 : V min < P[lp ] < Vedge < Per} (3.24) 

within which there are co-existing unstable / "wide " and stable / "narrow " solutions. 
There is also an open set in H 1 

{V>0 G H 1 : Vedge < V[M < V cr }, (3.25) 

where the only solitons are stable and "narrow". The terms wide and narrow refer, respec- 
tively, to solitons with frequencies in an interval near ( to the right of \x m in) or far ( to the left 
of (imin) the band edge, E* [54, 53] 

Finally, we state a 

Soliton excitation threshold conjecture (see also [64, 65, 33]V m in is an excitation thresh- 
old: If7*[ipo] < T^min, then ^(x, t) tends to zero as t — > oo ( LP, for some range of p > 2 
with the free Schrddinger decay-rate if ipo € H 1 and sufficiently localized in space) or in 
Lf oc , and scattering holds for ipo £ H 1 . 

Proof. Proof of Theorem 3.4. Part 1 follows from part 2 of Theorem 2.2, where we 
review results on the stability / instability of solitary waves. □ 

Proof. Proof of Theorem 3.5. This follows from an application of the sharp Gagliardo- 
Nirenberg inequality; see [61, 64]. Specifically, for any function / € i? 1 (R d ) we have 

i - ttS^V / iv/l 2 < I ( |v/| 2 - = H [f], 



\R\\ 



L' 1 



d 



(3.26) 
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where Ho denotes the conserved NLS/GP Hamiltonian for V = 0. Estimate (3.26) was used 
in [61] to establish, for V = 0, that if ip G H 1 and HV'olU 2 < \\R\\l», then NLS has a 
global in time _ff 1 bounded solution. It was further used in [64] to show that if, in addition 
we assume that \x\tpo 6 L 2 , then the solution decays to zero in LP, for range of p > 2 (and 
therefore in Lf oc ), □ 

Proof. To prove Theorem 3.5, note from (3.26) that 

1 - (j^) d ) / l y /l 2 ^ u \f\ J V I/I" < 3 - 27 ) 
Applying this inequality to a solution, i/j(x., t), of NLS/GP yields 

1 - (l^f)" ) / l V ^ x '*)! 2 ^ + J NAol 2 (3.28) 

For initial data, ipQ, in small iJ 1 neighborhood of a soliton with frequency near the band 
edge, we have H^oIIl 2 < ||-R||l 2 - Estimate (3.28) implies a uniform bound on \\ip(-, 
and therefore global existence (no singularity formation / no collapse). □ 

3.1. Finite gaps - results for focusing and defocusing nonlinearities. In this section 
we remark on extensions of our results to solitons with frequencies in finite gaps (gap soli- 
tons). For this, more general, discussion it is convenient to write NLS/GP and its nonlinear 
bound state equation in the form 

id t ip = -Aip + V(x)V> + g\4>\ 2a ip (3.29) 
( -A + V)u + gu 2a+1 =(j,u, (3.30) 

where we have introduced a parameter g to encode the (i) focusing / attractive (g = — 1) and 
the defocusing / repulsive (g = +1) cases. 

Focusing nonlinearity, g = — 1 : Our results of the previous section applied to solitons with 
frequencies in the semi-infinite gap, fj, < E*. The results on bifurcations of solutions from 
the spectral band edge can be extended to the case where E* is replaced by E e d ge , any band 
edge frequency. Here, we consider the case where the following two conditions hold 

1. The space of B - periodic solutions (— A+V)w(x) = E e d ge w(x) is one-dimensional, 
spanned by a function iy ec ; ffe (x), E e d ge is attained by the band dispersion function 
at k = 0. 5 

2. The inverse effective mass tensor, , is symmetric and positive definite. 

In this case, we have solitons centered about any point of symmetry of V(x), which in analogy 
to those described in Theorem 3.1, bifurcate from the left band-edge toward lower frequen- 
cies, into the spectral gap 

[1 1 ^ u(x, pi ) w (E edge - fX)^ io (x) F ( y/E edge - fj, (x - x ) J , 

E e d g e — \i > and sufficiently small. Here F satisfies the effective medium nonlinear 
Schrodinger equation (3.4), whose inverse effective mass tensor, A %1 is given by equation 
(3.5), with w replaced by w e d ge - Alternatively, this is (D 2 E n (ko))ij, the Hessian matrix of a 
Bloch dispersion function, E n , where EV^ko) = E e d ge , ko 6 8*. 

5 In dimensions d > 2 band edges may be attained at 5^ k £ B*; see [22]. In this case, the corresponding 
solutions are complex- valued and an extension of the present methods we use along the lines of [22] is necessary. 
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Defocusing nonlinearity, g = +1: Here, we consider the case where the following two 
conditions hold 

1. The space of B- periodic solutions (— A+V)iu(x) = E e d ge w(x.) is one-dimensional, 
spanned by a function w e d ge {x.), E e d ge is attained by the band dispersion function 
at k = 0. 

2. The inverse effective mass tensor, A 1 ^ = —B l \ is symmetric and negative definite. 
In this case, we have solitons centered about any point of symmetry of V(x), bifurcating from 
the right band-edge toward higher frequencies, into the spectral gap. 

Indeed, if we seek, along the lines of our previous analysis, soliton-like states with fre- 
quency: 

A* = Eedge — T(2 , 

our analysis near a band edge with negative definite effective mass tensor, — B 1 ^ , yields an 
effective medium soliton equation: 

d 

-Y, d v ^d V] F- lc uF 2 ° +1 =tF 

Thus, we can construct localized states for r < and p = E e d ge + |t|c 2 > E e d ge . 

Finally, we remark that all hypotheses concerning multiplicity of spectrum and curvature 
of band dispersion functions are verifiable in one space (d = 1) dimension. Thus we have 

THEOREM 3.6. Let V(x) denote a smooth, periodic and even potential. Consider any 
finite width, non-empty, spectral gap, — oo < a < b < oo, of —d 2 , + V(x). The band 
dispersion curvature at E — a is strictly negative and at E — b is strictly positive; see 
Appendix C. 

1. For focusing nonlinearity, g = — 1, centered about any point of symmetry of V, there 
exists a family of solitons of NLS-GP (3.29), which bifurcates from the zero solution 
with frequencies in the gap less than E = b. 

2. For defocusing nonlinearity, g = +1, centered about any point of symmetry of V, 
there exists a family of solitons of NLS-GP (3.29), which bifurcates from the zero 
state with frequencies bifurcating into the gap greater than E = a. 

These bifurcating branches have expansions and properties analogous to those described in 
Theorem 3.1 and Theorem 3.2. N.B. The results of this subsection indicate extensions to 
bifurcations into finite width gaps. In particular, for critical nonlinearities, we are able to 
analytically characterize the band-edge limit of the squared L 2 norm, V. Note however that 
the factor, £*, arising in finite gaps is associated with an excited Bloch state, i.e. a state 
Wed g e (x), which is not a positive ground state of the periodic boundary value problem. Since 
the estimate < 1, hinged on the result [34]: (m*) _1 < 1, which makes use of the ground 
state property (in particular positivity), we do not have an estimate on the size of in finite 
gap cases. 

4. Homogenization / multi-scale expansion . In this section we derive a formal mul- 
tiple scale expansion of solitons bifurcating from the band edge. In Section §5 we prove an 
error estimate, thus completing the proof of Theorem 3.1. 

Without loss of generality we choose coordinates with xq = 0. We seek a solution of the 
bound state equation (1.6), which bifurcates from the zero state at the band edge p, = E*, 
depending on a "fast" spatial scale x and a slow spatial scale 



y = € (x - x ) = ex, e < 1 
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(4.1) 



of the form 



He = E.« + e/xi + e 2 /i 2 + • ■ ■ (4.2a) 
u e (x) = e*C/ e (x,y) (4.2b) 

C/ e (x,y) = [/ (x,y)+ e (7 1 (x,y)+ e 2 C/ 2 (x,y) + ... . (4.2c) 

We also impose periodicity in x, i.e. 

C/ e (x + q J -, 2 /) = [/ e (x,y), j = l,...,d (4.3) 



Rewriting equation (1.6) by treating x and y as independent variables gives 

- (V x + eV y f U t + V(x)U t - e 2 C/ 2ff+1 = fi e U e 

Using the expansion (4.2) and the operator [see Eq. (2.6)], we obtain the following hier- 
archy of equations to C(e 4 ) 



O(e ) 


L*Uo 


= 0, 






(Die 1 ) 


LMi 


= (2V X 


■ V y +fXi)U 




0(e 2 ) 




= (2V X 


■ V y + /n)C/i 


+ (A y + M2 )C/ + [/ 2ff+1 


0(e 3 ) 


LJJ Z 


= (2V X 


■ v y + m)u 2 


+ (A y + /i 2 )C/i + (2ct + 1)C/ 2ct C/i + MaC/o, 


0(e 4 ) 




= (2V X 


■ Vy +fll)U 3 


+ (A y + /i 2 )£/ 2 






+ (2<7 


+ i)uS°u 2 + 


(2 ( j + l)(7C/ 2,T - 1 C/ 2 + /i4C/o 



where for each k > 5 we have: 

0(e fe ) : UV k = /j k Uo 

+ (2V X - V y + /ii)C/ fc _i(x,y) 

+ (A y + /i 2 )C/ fc _ 2 + J-fe [C/j (x, y), Mj : 1 < j < k - 2] (4.4) 

Note that is self-adjoint with a one-dimensional null-space spanned by w. In addition, fj, k 
is determined by a solvability condition of the form: 

/ifcM-)»^o(-,y)> + H-).^*(-.y)> = o, ( 4 -5) 

obtained by imposing orthogonality of id to the right hand side of (4.4). Here, F k denotes 
expression the sum of the last two lines on the right hand side of (4.4). Condition (4.5) ensures 
the existence of a solution to (4.4) which is periodic in x. 

We now implement this procedure at successive orders in e. In particular, we construct 
the terms Uj(x, y), < j < 4, as these are required in the proof of Theorem 3.1. 

4.1. Solution at each C(e fe ), k = 0,1,2,3,4. O(e ) terms: The O(e ) equation is 
solved by the choice 

Db(x,y) = «;(x)F(y), (4.6) 
where w is the periodic Bloch state associated with the band edge, as defined in Eq. (2.6). 

0(e 1 ) terms: The 0(e) equation for U\, by (4.6), becomes 

LMi = 2V x w • V y F + mw F (4.7) 
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Orthogonality of the right hand side of (4.7 to w implies fix — 0, from which we obtain 
(3.11): 

C/i = 2L; 1 [V X H ■ V y F(y) . (4.8) 



REMARK 4.1. To be completely systematic, we should add to the right hand side of (4.8) 
a term of the form Fn l (y)w(x), which is in the null space of ' L*, with -Fi^(y) to be deter- 
mined. Fih(y) is determined via the solvability condition for U3. Symmetry considerations 
lead to Fih(y) = (see the discussion of U3). We omit inclusion of this term to simplify 
the presentation. Note however that this degree of freedom is required at higher order. In 
particular, see the expression for t/ 2 (x, y) and the role of F 2 h{y) in the solving for Ui(x, y). 

C(e 2 ) terms: The C(e 2 ) equation for U 2 , by (4.6) and (3. 1 1) becomes 

LJJ* = tu(x) (A y + M2 )F(y) + ±d Xj d yj L- x [d Xi w){*) 9 yi F{y) 
+ W 2CT+1 (x)F 2ff+1 (y) 

= u.(x) A y F(y) + 4 V x • V y [ L' 1 [V x w] (x) • V y F ] 

+ /i 2 w(x)F(y) + w 2CT+1 (x)F 2,T+1 (y) (4.9) 

An equation for F(y) is obtained by imposing orthogonality of the right hand side of (4.9) to 
w(x). It is convenient to formulate the following 
PROPOSITION 4. 1 . Denote by the operator 

G(y)^A[G](x,y) 

= tu(x) A y G(y) + 4 V x • V y [ L^lV^-lfx) • V y G(y) ] (4.10) 

Then, 

(0,£*[G](-,y)) = d yi A^d Vj G(y) x (w,w). (4.11) 

Imposing orthogonality of the right hand side of (4.9) to w(x) and applying Proposition 
4.1 yields equation (3.4) for F = F(y, /j, 2 ): 

- d Vl A^d V] F{y^ 2 ) - lcif F 2 ° +1 (y,ii 2 ) = M2 F(y, M2 ). (4.12) 

Here, we consider only the positive decaying solution of (4.12), which by scaling and unique- 
ness can be expressed as 

F(y,fi 2 ) = IH^^XM'y;-!) (4.13) 
We can therefore, scale out |/i 2 | and henceforth assume /i 2 = — 1. 

Thus far, we have shown: 

To leading order, the slowly varying envelope function F(y) of the nonlinear bound state 
ofNLS/GP is comprised of a nonlinear bound state of the NLS equation (3.4) for a homoge- 
neous medium with effective mass tensor (A iJ ) [Eq. (3.5)] and effective nonlinearity 7 ff. 



17 



In Subsection §4.2 we show that F(y) is an appropriate scaling of R(y), the Townes soliton, 
the ground state associated with an isotropic homogeneous medium. 



We express the general solution of (4.9) in the form 

C/ 2 (x,y) = [/ 2 p(x,y) +EMx,y) = (7 2p (x,y) + w(x)F 2h (y), (4.14) 

where U 2p denotes a particular solution of (4.9) and w(x)F 2 / l (y) lies in the kernel of £* 
(recall L^w = 0), with -F^y) to be determined. 

Using equation (4. 12) to simplify the right hand side of (4.9) gives: 

C/ 2p (x,y)= J2 ^[(S^+Ad^L-'d^wWdy^Fiy) 

l<i,j<d 

+ w 2 ° +1 ( X )F 2 ° +1 (y) - w( X )F(y)] 

= L- 1 [( 5 l3 + W, /.. i), - A ij ) to(x) ] d Vi d yj F{y) 

+ L- 1 [ W 2ff+1 (x)- 7cffW (x)] F 2 ° +l {y) (4.15) 
or 

C/ 2p (x,y) = L-'xg^) d yi d yj F(y) + L^X^) F 2 ^(y), (4.16) 

l<i,j<d 

with X 2p j and X 2p ^ 2 given by the corresponding expressions in (4.15). To obtain (4.15), we 
use equation (4.12) for F(y), in terms of the effective mass tensor, (3.5) and effective cou- 
pling, (3.6). This is a consequence of the solvability (orthogonality) condition for equation 
for (4.9). 

C(e 3 ) terms: Using (4.14), we obtain the following equation for U 3 (x, y): 

LMs = 2V x -V y (C/ 2p +wF 2h ) 

+ ( A y - 1 + {2a + l)F 2a w 2a ) 2L- 1 d Xz w d yi F + fj, 3 w F. (4.17) 

Solvability of (4.17) requires orthogonality of the right hand side to w. Since all terms, except 
the last, on the right hand side of (4.17) are antisymmetric functions of x, we have fi 3 = 0. 
Thus, after substitution of the explicit expression for U 2p we have 

L*U 3 = 2d Xl w(*) d w F 2h (y) 

+ 2d Xl L: 1 Xg^)d Vl d yi d y] F{y) + 2d Xl L: l X 2pa { X )d yi F 2 °+\y) 
+ 2L- 1 d Xl w(x) (Ay -l)d yi F(y) 

+ 2w 2 °( X )L~ 1 d x Mx) d Vl F 2 ° +1 (y), (4.18) 
with summation over repeated indices implied. Thus, 

U 3 = 2L- 1 [\7 x w] ■ V y F 2h 
+ 2L- 1 

+ K 1 [V x w] • (A y - 1) VyF + 2w 2a L- 1 [V x w] ■ V y F 2a+1 ] 
= 2i; 1 [V x ™](x) . Vy F 2(l + U 3 (4.19) 
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VL^Kli] ■ V y d yi d yj F + VL^iX^-VyF 2 ^ 1 



0(e 4 ) terms: For U± we have 



L,Ua = (2V, ■ V„ + Mi) U s + (A y - 1)U 2 + (2a + 1)U 2 °U 2 
+ (j(2<j + l)^" 1 ^ 2 + fi 4 U 
= £*[F 2fe ](x,y) + (2<t + 1)™ 2ct+1 (x)F 2ct 
+ (A y - l)U 2p + (2a + l)U 2a U 2p + 2V X • V y U 3 
+ a(2a + l)Ul"- l Ul + fi 4 w(x.)F(y) 



(4.20) 



The operator £*[-](x, y), appearing in (4.20) is defined in Proposition. 4.1. Imposing or- 
thogonality of the right hand side of (4.20) and applying Proposition 4.1 gives the following 
equation for F 2 h'. 



LiF 2h (y) = (w^)- 1 [(w 1 (A y ~l)U 2p (- 1 y) + (2a+\)U 2a U 2p (- 1 y)) 
+ a(2a + l)(w,U 2a - 1 U 2 (-,y)) + 2(w, V x ■ V y C/ 3 (-, y)) 

= S(y), 

where L+ is the second order linear Schrodinger operator: 



L+ = -d Vi A"d Vj +1 - (2 ( T+l) 7cff F^+ 1 (y). 



^F(y) 
(4.21) 

(4.22) 



We now show that we can take fi 4 = 0. Equation (4.21) can be solved in L 2 (M. d ) for F 2 h(y) 
if and only if S(y) is L 2 — orthogonal to the kernel of L4. The kernel of has dimension d 
and is generated by translations, i.e. Kernel(L^) = span{9 w F(y), j = 1, . . . , d } [62, 39]. 
Since F(y) is even, the kernel of consists of functions which are antisymmetric in one 
coordinate direction. Moreover, it is easy to see that all terms in S(y) are symmetric and 
therefore orthogonal to the kernel of Lv:. Thus, we set /j, 4 = 0. 

4.2. F(y) is a scaled Townes soliton. Thus far, we have constructed the formal expan- 
sion (3.1) of (u e , n e ) through 0(e 2 ). The proof of its validity, in particular the error estimate 
(3.3), is given in Section §5. 

We conclude this section by relating the effective medium soliton F(y), which solves the 
NLS bound state equation with effective media parameters A 1 ^ and 7 e f f , to the unique ground 
state of the uniform-medium NLS equation, 



AR - R 



20-+1 



fxR, R>0, Re H 1 



(4.23) 



Let A = (A l; >) and A = diag(Xi, . . . , Ad) denote the diagonal matrix, whose diagonal entries 
are the eigenvalues of A. Let S denote an orthogonal matrix for which 



SAS 1 



A = diag(Xi, . . . , X c i) 



(4.24) 



Then, under the change of coordinates y z = A~?Sy, i*\(z) = F(y,— 1) solves 
Eq. (4.23) with fi = —1. By uniqueness up to translations, the solution to the isotropic 
NLS equation (4.23) is given by 



F(y) = Fi(z) = R(A-*Sy,-l) 



— ) 

Teff, 



! Sy,-l) 



(4.25) 



19 



Note that our expansion gives p, e = — e 2 + 0(e 5 ) (recall fi 3 = /i 4 = 0). As shown in the 
proof, we can in fact take 



Afi = fj, t = E, - e 2 . (4.26) 

Substitution of Eqs. (4.25) and (4.26) into the expansion (4.2), and using (4.6) and y = ex 
yields the leading order expansion of u(x, fi) displayed in Eq. (3.2). 

4.3. V [«(■, (j,)] near the band edge. To prove Theorem 3.2 and Corollary 3.3 we evalu- 
ate J \u e (x)\ 2 c?x, where u € is given by the two-scale expansion plus error term (3.2), of The- 
orem 3.1. We obtain (again recalling the choice of coordinates, so that y = e (x — xo) = ex) 

J |u £ (x)| 2 dx = ei J |[/ (x,ex)| 2 +2eC/ (x,ex)[/i(x,ex) 

ao ai 

+ e 2 (|C/!(x,ex)| 2 +2[/ (x,ex)[/ 2 (x,ex)) dx + C(e 3 ) 

V v ' 

Q 2 

= d [T + X\ + 11 + G(e 3 )] (4.27) 

Each of the three terms on the right-hand side will be treated below using the following 
general averaging method 

LEMMA 4.2. Let p(x) be periodic on the lattice V having the fundamental period-cell 
B. Let B* denote the dual fundamental cell (first Brillouin zone) which spans the dual lattice 
r*. Assume that X^ker* l-Pk| < °°> where {pu} denotes the set of Fourier coefficients of p. 
Let G e L 1 (R d ) n G°°(R rf ). Then, as e -> 



/ p(x) G(ex)dx - / p( X )dx ■ f G{y)dy = G(e°°) x V \ Pw \ = G(e°°) 
Js. d Jb Jm d ker , 



where the cell average § B is defined by £ B p = jj^ J B p. 

Proof. Proof of Lemma 4.2. p(x) has the Fourier representation 



ker 

Therefore, 



p(x) = V] p k e ik x , where p k = f e lk x p(x)dx . 



(4.28) 



t d [ p(x) G(ex)dx = e d V p k / e* kx G(ex)dx 

= p / G(y)dy + V Pk [ e^G(y)dy 



: / p(x)dx / G(y)dy + V Pk G 
Jb Jr* 0/ker , 



k 

2?re 



By smoothness of G, for all q > 1 and £ £ M. d there is a positive constant, ra,q, such 
that |G(£)| < ?'G.g(l + The required estimate of the remainder term follows. This 

completes the proof of the Lemma. □ 
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We now proceed with proof of the Corollary 3.3 by evaluating the terms 2|, j = 0, 1, 2 in 
Eq. (4.27). 

Claim 1: 

T = e- d C* T cr + 0{e°°), where (4.29) 



Proof. By Eq. (4.6) and Lemma 4.2 one has 

2g= [a dx= [ w 2 (x) F 2 (ex) dx = e~ d / iu 2 (x) cix f F 2 (y) dy + 0(e°°). 

J JR d JB JR d 

Using expression (3.8) for F(y) as a scaling of R(y; — 1) we get 

/ F 2 (y) dy = (7ef f )-- / i? 2 (y;-l)dy 



2 

m* 2 Vcr , 



2cr+2 



□ 

Claim 2: If = 0(e°°). 

Proof. We proceed similarly by using (3.11) and Lemma 4.2. We obtain 

X\ = [ati dx = 2e V / w(x) F(ex) • 2L; 1 (x)d % F(ex) dx 

= 4e 1 - d £ ^(x)!- 1 (0 aj tc)(x) dx J F{y)d Vj F{y) dy + e 0(e°°) = 0(e~) , 



since f Rd F(y)d vt F(y)dy = 0.0 
Finally, we turn to 2J. 



Claim 3: 



2 2 = e" d - e 2 



4 E / l £ * 1 t^^W] I 2 dx / I V(y) I 2 rf y 

£w 2 (x)dx J d n F(y)S(y) dy 



0(e°°). 

(4.31) 

where S(y) is explicitly displayed in (4.21). 
Proof. 

T 2 = Ja 2 dx = e 2 J [|(7i(x,ex)| 2 + 2[/ (x, ex) (U 2p (x, ex) + w(x)F 2fe (ex) ) ] dx 

= 22 iQ + 22,6 + 1-2.C ■ 
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and, by Lemma 4.2, X\ a is given by 

T-la = e 2 / |C/i(x,ex)| 2 dx = 4 e 2 J \L? [d Xj w{x)} d y .F{ex) | 2 dx 

= ^ d { 4 E / l^ 1 K^W] I 2 dx ■ l I V(y) \ 2 d y + 

Concerning l 2 ,b, we assert the following: 

x lb = 2e2 / w{x)F(ex)U 2p {x 7 ex) dx = 0(e°°). 

□ 

Proof. To prove (4.32) we note that U 2 [Eq. (4.16)] is of a sum of terms that have the 
factored form 

U 2p {x, ex) = ]T G 3 (ex) ■ L^P^g^x) . (4.32) 

i 

Here, P- 1 denote the projection onto the orthogonal complement of w in L 2 eriodic (B); see 
(2.8). Substitution of (4.32) gives 

Zi,b = 2e 2 V I w(x) L~ 1 P ± g j (x) ■ F(ex)G j (ex) dx 
which by Lemma 4.2 implies 

Z 2, 6 = Y,£ w W L * lp± 3iW dx j F(y)Gj(y)dy + 0(e°°). 



Since P- 1 commutes with functions of L* and P^w = 0we have 

(w(x.),L- 1 P ± g j ) = (w^-'P 1 - P ± g J ) = (P^L; 1 P ± g J ) = 0. 
It remains to calculate I 2)C . 

l 2 , c = 2e 2 J U {x,ex)U 2h (x,ex) dx 

= 2e 2 J w(x)F(ex) ■ w(x)F 2h (ex) dx 

= 2e 2 - d £w 2 -J F(y) F 2h {y) dy + 0(e°°), 

= 2e 2 - d £w 2 - J F(y) (Lf)- 1 S(y)dy + 0(e°°) 

2e2 ~ d f B u ' 2 ' I ( L +y lF ^ s ^ d y + 

-F(y) + y ■ V y F(y)) S(y) dy + 0(e°°) (4.33) 



IB 

-2 
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Here we have used the relation 



[L$) F=d ll2 F{-,n 2 ) 



1 fl 



ft2=— 1 



-F(y)+yV y F(y] 



which follows from differentiation of the equation for F = F(x; ^2) with respect to see 
(4.12) and (4.13). 

Therefore, summing up the terms we have 

k(x)| 2 dx = ei [XI + X{ + X\ + 0(e 3 ) ] 



o 2 



L 

d 



4 E / l^ 1 [o*M*)] I 2 dx / I d vi F (y) I 2 d y 

]= lJB JR d 

dx J (^F(y) + y • V y F(y)) S(y) dy 



+ e>(e°°) 



= (e 2 )^"' (C*Pcr +e 2 Ci* +0(e°°)) 

Recall 5(y) is displayed in (4.21). □ 

This concludes the proof of Theorem 3.2. 

5. The error estimate (3.3) and conclusion of the proof of Theorem 3.1. In this sec- 
tion we prove Theorem 3.1. For ease of presentation, we focus on the cubic (cr = 1) one- 
dimensional case (d = 1): 

( -dl + V(x) ) u - u 3 = fj, u 

The proof carries over to the more setting in the statement of Theorem 3.1. After the proof, 
we indicate the modifications required for the proof to go through in general dimension d = 
1,2,3; see Remark 5.2 below. 

We shall construct a solution (u, E) = (u e , fi e ), using the formal multiple scale expan- 
sion of Section §4 : 



u t = eU e {x) 



^ = E, -e 2 



L k=0 



(5.1) 



(5.2) 



The expansion includes an error term, e 3 C/|(x), which must be estimated. The equation for 



e 3 UI is: 



( -dl + V(x) - 3e 2 C/ 2 (a;, ex) - E*+e 2 ) Ul{x) 



,2 ne 



R e [Uj;0<j< 4,U£(x)] 
e 2 R e o [U r ,0 < j < 4] + e 3 R\ [U r ,0 < j < 4] U£ 
+ e 5 i? 2 [[/,;0<j<4]([/|) 2 + e s ( C/|) 3 , 



(5.3) 



where R € k [Uj ; < j < 4] denotes the coefficient of the k power of t/| , and is a polynomial 
in the previously constructed functions Uj, j = 0, 1, 2, 3, 4. 
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The scaling of the error term in (5.1) is motivated as follows. Formally, the correction to 
the leading order sum in (5.1) will be of order e 5 . In our analysis, we find that the frequency 
components of the corrector (to the truncated multiple scale expansion near the band edge 
are of order 0(e 3 ) . Therefore anticipate this result in (5.2). We will in fact show that for 
s > d/2, \\U§\\h s is bounded uniformly in e. This implies the error bound of Theorem 3.1. 

In particular, 

Rl = 2d x d y U 4 + (d 2 - 1)U 3 + 3U 2 + 6[/ot/iE/ 2 + Uf + 0(e) . 

Our goal is to estimate C/| and to do this we employ the spectral (Floquet - Bloch) decompo- 
sition of the operator — A + V(x). 

5.1. Floquet-Bloch Theory and the Bloch transform. See the references [23, 47, 36] 
for basic results on the spectral theory of operators with periodic coefficients. 

Assume V(x + 2n) = V(x). For each k e T = [— i, i] we seek solutions of the 
eigenvalue equation for the operator (— d 2 + V(x)) of the form: 

u(x; k) = e lkx p(x; k), p(x + 2tt; k) = p(x; k), x e K 

This yields the periodic elliptic eigenvalue problem for p(x; k) 

(-(d x + ik) 2 + V{x)) p(x;k) = Ep(x;k) 7 p(x + 2n; k) = p(x; k) 

For each k E T the spectrum is discrete give rise to eigenpairs (E m (k) , p m (x; k)) m >i, and a 
complete orthonormal set {p m (x; k)} in L 2 er with respect to the inner product: 

(fi9)Ll er = ^ f{x)g{x) dx. 

(E n (k) 7 u n (x;k)) 7 n>l, fee [—1/2, 1/2] are solutions of the eigenvalue problem: 

(-dl + V(x))u n (x;k) = E„(k) u n (x; k) 
u n (x + 2tt; k) = e 2 ^ lk u n (x] k), i£R, 

where k i-> E m (k) sweeps out the m*' 1 spectral band, and yield a complete set of states in 
L 2 (]R);see (2.5). 

N.B In this section we assume that w(x) is normalized, (w, w) = 1. Thus, w(x) is the unique 
normalized ground state of the periodic boundary value problem and 

(pi(z;0),£i(0)) = (p 1 (x;0),E 1 (0)) = (w(x),E t ) 

Furthermore, for each k 6 T, the set {p n (x; k)} is an orthonormal set in L 2 er ( [0, 2tt) ). 
Introduce the Gelfand - Bloch transform (T<fi)(x; k) = 4>(x; k), and its inverse T~ x : 

(T0)(x; k) = j)(x; k) = J2 e im " x 4>(k + m) (5.4) 

mGZ d 

(T~ 1 0)(x) = / e lk - x 4>(x;k)dk, 

where 4>(k) denotes the Fourier transform of 4>(x). Clearly we have 

4>(x + 2ir; k) = (j>(x; k), and 4>{x\ k + 1) = e~ lx 4>(x; k) 
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One can check that 



TT~ l = Identity on L 2 (R). 

Another important property of T is that it commutes with multiplication by a periodic 
function: 

/(i + 27r) = /(a:) (Tfg)(x;k) = f(x)(Tg)(x) 

Since <fi(x; k) is 2tt— periodic in x, we have 

4>(x;k) = ^{p m (-;k),4>(-;k))p m (x;k) (5.5) 

m>l 

We conclude this subsection with some basic definitions and results required below; see, for 
example, [22] and references cited therein. 
Theorem 5.1. 

1. There exist positive constants Ci,C2 band dispersion functions E n (k), n > 1 satisfy 
the bounds 6 



c\ n 2 < E n {k) < c 2 n 2 , \k\ < 1/2 



(5.6) 



see [20, 31 ]. 
2. The mapping 



0) >-> 



</>n(fc) 



;s an isomorphism of H S (M 1 ) with X s = L 2 (T X ; l 2,s ), with norm: 

/ &]T(i + M 2 )" 





2 






\ / n>l 





.Y- 



(5.7) 



J. Moreover, there exist positive constants C\ , C2, swc/i f/iaf we have the norm equiva- 
lence 



< 



,k) , Pn(-,k) 



: C 2 



Lffs 



(5.8) 



.v- 



4. Assume 4>,ip € H s (R d ). 

(i) Ifs > q+d/2, then cj> € C^(R d ), the space of C q functions, f, with \d a f{x) \ -> 
as x — > 00, |a| < q. 

(ii) Ifs > d/2 then H s is an algebra, i.e. <fiip £ H s and \\ <fiip \\h' < C II 011 if s IIV'II-H 3 - 
REMARK 5.1. The bounds ( 5.6) are well known; see [20, 31 ]. To prove the isomorphism, 

recall the operator = —A + V — E* > 0; see (2.7 ). Standard elliptic theory implies 



6 In dimension d, n 2 is replaced by n d . 
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\ 2 H * ~ + l*)H\\1 



Li 4>\\l^ defines a norm equivalent to the H s norm. Furthermore, by (5.5) 

I e* Y Mk) (1 + E* - E n {k)) s Pj (; k) 
\Mk)\ 2 \1 + E, -E n (k)\ s dk 
|0„(fc)| 2 (l + |n| 2 r dk 

2 

X s 



L- 



E 

n>\ 

E 



n>l " 



n>l " 



<f>n(k) 



(5.9) 



5.2. Corrector equation and localization in Bloch variables. In this subsection we 
express the equation for the corrector 



^(x) = Ul(x) 



(5.10) 



in Floquet-Bloch variables and, in particular, decompose this equation into spectral compo- 
nents near and away from the band edge E* ; see, for example, [ 14, 2 1 , 22] . 

Applying the Bloch transform, T, to (5.3) we obtain an equation for $? e (x; k) = (T^ e )(x; k): 

- (d x + ikf + V{x) -E,+ e 2 ] (T* e )(x; k) 
- 3e 2 w 2 (x)T [F 2 (e-)<T(-)] (x;k) 

= e 2 (TR e )(x;k), (5.11) 

where e 2 R e is defined in (5.3). Here, we have used that Uo(x, y) = w(x)F(y). 
Now \P c (x; k) is periodic in x. Therefore, 

oo 

&(x;k) = & m (k)Pm(x;k), & m (k) = ( Pm (-;k), * e (-;fc)> 

m— 1 

We introduce a decomposition of \& e into spectral components near the band edge E\ (k = 
0) = E* {low frequencies) and spectral components away from E* ( high frequencies) as 
follows. Let 1^ denote the characteristic function for the set A and define 

X(a < k < b) = l{ fe:Q < fc < 6 } . 

Express fy e (x, k) as 

&(x;k)= X (\k\ <e r )^t(k) Pl (x;k) 



\P? (x;fc) 

X {e r <\k\<2- 1 )^ 1 (k) Pl (x;k) + £ % n (k) p m (x; k), 

m>2 



(5.12) 



where r is chosen to satisfy 



< r < 1 . 
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(5.13) 



Using the inverse Bloch transform, we obtain 

high 



Taking the inner product of (5.1 1) with Pj(-; k), we obtain: 

[E j (k)-E, + <?] &j(k) - 3e 2 (p,(-;fc), w 2 (-)T [F 2 (e-)^(-;k)} (■;*)) 
= e 2 ( Pj (-,fc) , (Ti?*)(-,*)) = £ 2 (Ti^(fc), J > 1 

(5.14) 

The system (5.14) can be viewed as two coupled systems for the low and high frequencies: 

*W*) = X(I*I <0*!(*) and 

HighW = (x(e r < 1*1 < 2- 1 ) * e (fc) , {f^)b>2 ) 

(5.15) 

Low frequency components: 

[^(ftJ-B.+e 2 ] %, iow (*) 

- 3c 2 x(|A|<c r ) (Pi(-,k), w 2 (-)T[F 2 (e-)n,i ow (-)] (■;*)> 
= 3e 2 X (|fc|<0 ( Pl (-;fc), ^ 2 (-)T [F 2 (e)*^(-)] (s*)) + ^ 

(5.16) 

High frequency components: 

[ E^k) -E*+e 2 ] X (e r < \k\ < 2" 1 ) 

= 36 2 X ( £ r < |*| < 2- 1 ) ( Pl (.;fc) , u; 2 (-)T [F 2 ( £ .)* e (0] (■)) + ^ lWfl /. 

(5.17) 

[£,-(*) -£*+e 2 ] §|(fc) 

= 36 2 (^(.;fc), W 2 (-)T[F 2 ( £ .)* e (-)](-;fc)) + e 2 Rl high , j>2 

(5.18) 

Here, R\, ow - 





(Rj,high)i>i are given by 






= x(|fc|<e r )(pi(-,fc), (Ti? e )(-,*)} 


(5.19) 


Rl,high 


= X(e r <|fc|<2- 1 )(p 1 (-,fc), (Ti? £ )(-,fc)) 


(5.20) 


^j.high 


= fe(-,fc), (Ti? e )(-,fc)), i>2, 


(5.21) 



where i? £ is defined in (5.3). We study the system for ^\ i ow (k), ^\ igh {k) usm g the fol- 
lowing 

Lyapunov-Schmidt reduction strategy: 

1. Using the implicit function theorem, solve the infinite system of high frequency 



component equations for ^% lgh as a functional of %\ <low : ^ e high = V e high 
with an appropriate bound on this mapping. 

27 



low 



2. Substitute &% igh = ^high ^1 low mt0 (5-16) to obtain a closed equation for the 
low frequency components, which is solved via fixed point iteration. 
We now embark on implementing this strategy. Our first step is to rewrite the low frequency 
equation (5.16) in appropriately rescaled variables. 

5.3. Closure and rescaled of low frequency equation for ^\ low . With a view toward 
obtaining a closed equation for , ow , we begin with several observations. 

1. From our formal multi-scale construction, we expect ^ , on] (i) ~ $(ex) w(x). This 
motivates the following 

Ansatz: Seek the low frequency components in the form 

*U(*0 = X(\k\<e r ) ^*(*), (5.22) 

Thus, 

m ow {x;k) = X (\k\<e r ) 7 6 (7) P^ k )- 
Using the definition of T -1 and that 

p 1 (x;k)=pi(x;0) + O(k) = w{x)+0{e r ) 1 \k\ < e r we have: 

= Hex) w(x) + 0(e r ). 

2. Note that for |fc| < e r , 

E^k) — E* - i^S 1 (0)fc 2 - ife(ft) fc 3 , < k < e r 

Thus, 

(E x (k) — E* — e 2 ) % low (k) = Qfe(o)fc 2 - 6 2 ) X (\k\ < old (?) 



(5.23) 



This and the Ansatz (5.22) suggest the scaling 



r 

In this scaling (5.23) becomes 



k 

k = - (5.24) 



= e 2 ( ld 2 E 1 (0) K 2 - 1 ) x (|«| < e^ 1 )-^) 



[ e 3r X (|«| < e r - 1 )i$( K ) 
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3. Consider the last term on the left hand side of (5. 16). We have 

- 36 2 X (|fc|<0 ( Pl (.;k),w 2 (-)T[F 2 (e-)n,io W (-)] (s*)) 
= - 3e 2 X (\k\ < e r ) ( Pl (-, k) , w 2 (-)T [F 2 (e) x (|V X | < O *(e) «;(•)] (■; k)) + 0{t 
= -3e 2 X (\M < O (px(-; fc) , w 3 (-)T [F 2 (e-) X (|V X | < O *>(*)] (•; fc)) + 0(e r+2 ) 

= -3e 2 | W 4 ^- X (|fc|<e r ) [F 2 (e)x(|V 3; |<eO$(e)]"(fc)+ 0(e r+2 ) 
= -3e 2 JwUx-x (\k\ < e r ) \ [F 2 X (|V„| < e^ 1 ) + 0(e r+2 ) 

= - 3e 2 /V • x (N < e'- 1 ) - [F 2 X (|V„| < e''" 1 ) $]" («) + Q(e r + 2 ) 



r+2\ 



Use of (5.22) and (5.24) in (5.16) yields the following: 
Closed / Rescaled low frequency equation for $(«): 



x(N<c r 



7? 6 , , , ( K- $ ■ ■ 

-"rescaled low V * high 



(5.26) 



Here, £i : [—1/2, 1/2] — >• [.E*, £1(5)] , fe -Ei(fc) denotes the band dispersion function 
for the first spectral band and 7 e ff is given by the expression in (3.6). As we have assumed 
Pi(x, 0) = w(x) to be normalized, 7 c ff = J 2 * w 4 . 

We summarize the arguments of this subsection, which lead to the system we'll study: 

PROPOSITION 5.2. The coupled system consisting equation (5.26) for the rescaled low 
frequency components, < e r_1 (|fe| < e r ), and the high frequency equations, 

(5.17) and (5.18) is equivalent to the original system. 

5.4. Proof by Lyapunov-Schmidt reduction. We estimate the right hand sides of the 
high frequency equations (5.17) and (5.18). 

PROPOSITION 5.3. Let s > d/2. For some positive constants C\ and C2 we have 



3e 2 ( Pj (-;*), w 2 (-)T[F 2 (e-)^(-)} (■;*)) \\ xs < d e 2 ||* e |U 



R 



j,high 



< a 



X* 



(0(e°°) 



e 3 ||* e| 



<*Wtx» 



e 8 Wf x , 



PROPOSITION 5.4. The system for ^ e high can be solved in terms of 
^1 low = X{\k\ < O^i low an dwe have the estimate: 



high 



low 



< C (<9(e°°) + e 3 - 2r \\ X (\k\ < e r M,iow IU» + e 5 ~ 2r \\ X (\k\ < e r M, low \\ 
+ e 8 ~ 2r \\x(\k\ < e r M ,i ow fx 



(5.27) 
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Proof. Consider the system for ^% igh , (5.17-5.18). The result follows from direct esti- 
mation using: 

| Ej(k)-E* +( 2 | > c> 0, j > 2, 

| E^k) -E*+e 2 \> e 2r , e r < \k\ < 1/2 (5.28) 

and applying the implicit function theorem. □ 

Use of the Proposition 5.4 in the rescaled low frequency equation we obtain 
Proposition 5.5. 

i^'(0) K 2 + l) x(\K\<* r - l )*(*) - 3 7eff [F 2 x (|V y |<6 r - 1 )$]» = H' 
where H e = x (\ K \ ^ e r H e and satisfies the bound 
II < U&ym + O (\\d 3 E 1 \\ 00 e 3r \\ X (|«| < e^ 1 ) $|U=,« ) 

+ e CT (|| X (M < e^ 1 ) $|| ia .. + || X (|«| < e^ 1 ) ) (5.29) 

w/zere ct > ant/ C/ sym € H^,m- 
Equivalently, we have 

( -9, Ad y - 3 7off F 2 + 1 ) X (|V„| < e^ 1 ) $ = X (|V„| < e^ 1 ) 



(5.30) 



w/zere ^4 = hE"(Q) and 



2 

II # £ \\ H s < \\Q sym \\ H s + O (H^IU e 3r Hx (|V„| < e^ 1 ) ) 
+ ^ (II X (|V V | < e^ 1 ) + || X (|V V | < e- 1 ) *|&. ) 

We now complete the proof. Denote by the operator 

Li = - d y Ad y - 3 7cff F 2 (y) + 1 (5.31) 

and 

Xe = X (|V y | < e 1 - 1 ) , x7 = 1 - X, = X (|V„| > e^ 1 ) (5.32) 
We recall (see (5.13)) 

< r < 1 . 

Equation (5.30) for $ can be written as 

Xei+Xe $ = (5.33) 

Since F is chosen to be centered at local extremum of the symmetric potential, V(x), we 
have that the mapping 

$ ^ X e 
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maps H* ven to itself. 

We claim that for some e , if e < e , then the operator \ t Xe ■ Hsym ^sym 
has an inverse with norm bound which depends only on eo. Thus, for < e < ei, we can 
reformulate (5.33) as 

and show by fixed point iteration that for some < t\ < eo, sufficiently small, equation 
(5.33) has a unique H s+2 solution, which is bounded uniformly for e < e\. This then implies 
Theorem 3.1 

Therefore, the proof boils down to establishing the invertibility of Xe Xt '■ H-lym ~~ ^ 
-Hsym- We first prove that : — > H^ ym has a bounded inverse. 

Now the the operator L+, acting in L 2 (R 1 ) has one-dimensional kernel, spanned by the 
function d y F. To see this, differentiate the equation for F(y): 

- d y Ad y F + F - 7cff F 3 = 

and obtain L+F' = 0. Moreover, Ker(L+) = span{i 7 "}, since the eigenvalues of a Sturm- 
Liouville operator are simple. Since F' is odd, is an invertible and bounded map from 

Finally, we turn to the invertibility of Xe Xc '■ Hsym ^sym f° r e sufficiently 
small. We begin by expressing Xe Xe as a perturbation of L+: 

X, L% Xe = L% + Q € 

Q e = -(TeL* + L^— c ) + x^L^Te 

Therefore it suffices to prove that 

Ll + Q, = L% (I + {L%Y X Qe) 

has a bounded inverse defined on H^ ym . A bounded inverse 

(Li + Qty 1 = (i + {L+y 1 q e y 1 (L^y 1 

exists provided the norm of Q e can be made smaller than one, by choosing e suffi- 

ciently small. 

(Lir'Q, = -(L^T^Lix, XT 

Concerning the second term, the mapping / i-> x7f maps H s to H". If s > 0, the operator 
norm tends to zero as e — > 0, by explicit calculation using the Fourier transform. 

Finally, consider the mapping / h-> xl We prove that this mapping is 

bounded from H k to H k ~ s , for any S > 0. We see this as follows. Denote by (y) = 
(1 + |y| 2 ) 5 and therefore the operator (D) a is defined by 

(D) a f = J e*-* (k) a f(k) dk 

Now, for any a > 0, we write 

(L^TeLi = (L^)" 1 (D) a ■ {D)- a x e ■ Li 
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and estimate the norm as follows: 



< H(if)- 1 (D) a \\ H .-^ H .-, \\(D)-*Xe\\H>-^H>-' ■ \\L$\\ B .-^ H . 

Note that the first and third factors are bounded independently of e. We claim that J| (D)~ a x e \\H s 
as e — > 0. To see this, calculate as follows: 

\\(D)- a x t f\\h = j^r 2a i {lKl <^ } ^) T \m\ 2 d K 

< e 2a(1 - r) !l/|| 2 ff . 

Thus for any a > 0, we have || (D)~ a XeWn^^H^ — > This completes the proof of Theorem 
3.1. 

Remark 5.2. General spatial dimensions d > 1: The proof given readily extends to 
general dimension d > 1. One works in spaces H s (R. d ), s > d/2. X s is constructed taking 
into account the behavior of the dispersion functions, E n (k), in dimension d. As before, 
we choose F(y) = F(e(x — Xo)), with Xo a point of symmetry ofV(ic). The kernel of 

= — d Vi A' l: 'dy j + F — 3F 2 (y) has dimension d and is generated by translations, i.e. 
Kernel(L^) = spa.n{d Vj F(y), j = 1, . . . , d } [62, 39]. Since Kcrncl(L^) is orthogonal to 
H% ven {R d ), L\ is invertible mapping from H s+2 (R n ) to H s . 

6. Numerical computations in the semi-infinite gap. Our analytical results apply to 
solitons with frequencies in a spectral gap, which are also sufficiently close to a spectral 
band edge. In this section we present the results of numerical computations corroborating 
the rigorous asymptotic results near the spectral band edge, but also illustrating their approx- 
imate validity further away from the band edge, well into the spectral gap. The details of the 
numerical methods are discussed in Section §6.3. 

The particular rigorous asymptotic results we explore numerically in detail are: 

1 . the asymptotic structure of soliton's lying near the edge of the spectral gap (Theorem 
3.1): 

u(x, n) w {E* - fx)' w(x) F (e(x - x )) , 

where xo is a local extremum of V(x). 

2. The asymptotic behavior of the soliton (nonlinear bound state ) power, V(p), along 
minima- and maxima-centered solitons as p approaches E*; see Theorem 3.2 and 
Corollary 3.3. In particular, in the critical case a = 2/d we have 

V(p) w C* V cr (6.1) 

We focus on the one-dimensional NLS/GP equation (1.1) with critical nonlinearity and 
periodic potential governing tp(x, t) and nonlinear bound states: ip(x, t) = e~ lA1 *u(x, p). 

d= l,cr = 2: 

id t i> = -d 2 J) + Vo cos 2 (2ttx)iP - (6.2) 
pu = —d 2 u + Vq cos 2 (2irx)u — |u| 4 u (6.3) 

Vo is the variation or contrast of the potential. 
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We have observed similar results to those presented below for the two-dimensional criti- 
cal NLS/GP with periodic potential: d = 2, cr = 1: 

id t ij; t = - (d 2 x + dfj i> + y [cos 2 (27rx) + cos 2 (27ry)] V - \ip\ 2 ip ■ (6.4) 

liu = - (d 2 x + <9 2 ) u + y [ cos 2 (2tt:z;) + cos 2 (27ry) ] u - \u\ 2 u (6.5) 

Theorems 3. 1 and 3.2, and Corollary 3.3 apply to (6.3) with states centered at a minimum: 
.To = .25 or maximum: xq = 0. These results also apply to (6.5) with states centered at a 
minimum: xo = (.25, .25), at a maximum: Xo = (0, 0), and at a saddle point: xo = (.25, 0) 
or x = (0, .25). 

6.1. Soliton profiles: asymptotic theory and computation. Figures 6.1 and 6.2 dis- 
play nonlinear bound state profiles of the one-dimensional NLS/GP equation (1.1) for values 
of /i in the semi-infinite gap of the Schrodinger operator: <9 2 + Vq cos(Kx), i.e. 

[j,€ (-oc,E*), ii < E* = E*(V ,K) 

both near and far from the band edge. 

Plots Al, A2 and A3 in Figure 6. 1 display the case of soli tons, centered at local minima of 
the potential with, from left to right, frequency /i approaching E« , at distances E«—fi= 10, 1 
and 0.01, respectively. Plots Bl, B2, and B3 in Figure 6.1 correspond to the case of solitons 
centered at local maxima of the potential. 

We first note that the figures show the expected trend toward increased localization as 
fi < is decreased. For /i large and negative the solitons centered at maximum or minima 
approach a scaled V = soliton; see (1.13) and 

Our main analytical results apply to solitons whose frequencies lie near the band edge, 
although numerical studies indicate their approximate validity some distance away from the 
band edge. 

Theorem 3. 1 implies that nonlinear bound states are, to leading in order in the distance to 
the spectral band edge, a product of a linear Bloch state with band edge energy and a soliton 
in an effective homogeneous medium; see (3.2) and (1.13): 

u(x,fj,) ps w(x) F(y) 
( E — \ 

F(y) = ( secbi ( 2y/m,(E, - //) (x - x ) ) . (6.6) 

The centering of the soliton is xq, a point of symmetry of the potential, V(x). 

The maximum of the Bloch modes w(x) (normalized to be positive and with unit mass) 
occurs at the minimum of the potential. However, depending on the centering point Xq, 
F(y) has a maximum (minimum) at the potential maximum (minimum). Thus, bound states 
centered on a potential minimum are approximately a product of functions that peak at the 
same values of x yielding a more peaked bound state; compare the top and bottom panels of 
Fig. 6.1 within - fi = 10. 

Figure 6.1 (A2,B2) shows that for E* — fj, = 1 the bound states have discernible oscil- 
lations about a positive envelope, reflecting the solutions leading order behavior (6.6). These 
oscillations can be understood as a result of the "underlying" Bloch modes in Eq. (6.6). Here 
as above, the asymptotic theory appears to capture the structure of the bound states even when 
/i is not very close to the band edge. 
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We note as well, for the soliton centered at the potential's local maximum, a transition in 
the profile from single-humped to double-humped (having a dimple at x = 0) as ft decreases 
through /i = fi#, the value at which V[u(-, fi)], along the branch of solitons, centered at a 
local maximum ofV, achieves its maximum; see figure 6.3. A related observation is made in 
[1]. 

Comparing Figs. 6.1 (A3) and 6.1 (B3) shows that near the band edge (E* — (J. = 0.01) 
there is almost no visible difference between the bound states centered at potential minima 
and those centered at potential maxima. This is clear from Eq. (6.6), since in this regime 
F(y) ~ decays only on a length scale much larger than the period of V(x) and thus, for both 
maxima- and minima- centered solitons, u(x) ~ constant xw(x). 

Figure 6.2 shows a direct comparison between the asymptotic theory, i.e., the leading 
order solution Eq. (3.2) near the band edge, and the "actual" bound state profiles, computed 
by solving the bound state differential equation (6.3) with high accuracy. 



A1 A2 A3 




FIG. 6.1. Leading order asymptotic profiles obtained via Theorem 3.1 (blue, solid) centered on the potential 
minimum (top panel) and potential maximum (bottom panel) for (Al, Bl) E* — p = 10 , (A2, B2) E* — p, = 1, and 
(A3, B3) E f — p = 0.01. Also shown are the scaled (for plotting purposes) potential V(x)/Vo (red, solid). Block 
wave w(x) (green, dashes), and rescaled homogeneous ground state F(y) = F (e(x — xq)) [Eq. (6.6), black, 
dots)]. For clarity only a small portion of the domain is shown and the x-axes are zoomed in as E* — p decreases. 
Geometric shapes correspond to those depicted in Figs. 6.2 and 6.3. 

6.2. Effective mass and the power curve \i >->• V[u{- , fj,)]. Each plot in Figure 6.3 
shows two curves of bound-state power, "P[it(-, (i)], in the semi-infinite gap as a function of 
H for (6.3) with Vb = 10, K = 2n. The solid (blue) curve corresponds to the variation of 
V\u{-, /j.)] for the family of solitons centered at the potential's local minimum and the dashed 
(red) curve for the family centered at the potential's local maximum. 

We observe the following: 
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XXX 

FIG. 6.2. Bound-state profiles computed using the Renormalization method [Eq. (1.6), blue, solid] compared 
with the leading order asymptotic theory [Eq. (3.2), red, dashes] for the same parameters as in Fig. 6.1. 

1 . Panel ( A 1 ) of Figure 6.3 shows the variation of V over a wide range of fx in the semi- 
infinite gap. As fx —> — oo, i.e. in the semi-classical limit the power approaches the 
homogeneous (V = 0) power of the ground state, V cr = -^7r ss 2.72 . 

2. Panels (A2) and (A3) of Figure 6.3 show, as predicted by Corollary 3.3, that as \x — > 
E* the ■p[u(-,/i)] approaches the value (*V C r ~ 2.2 strictly less than Vcr ~ 2.72. 
Here, ^ as 2.2/2.72 as 0.8. This is true for solitons centered at either potential 
minima or potential maxima; see [Fig. 6.3 (A3)]. 

3. Although the asymptotic behavior of the V for maxima and minima- centered bound 
states is the same, across most of the gap bound states centered on lattice minima 
(resp. maxima) have power below (respectively, above) V C r ~ 2.72. 

4. Panels (Al) and (A2) of Figures 6.3 show a transition in the slope of the minima- 
and maxima- centered power curves at the same value [i = As discussed 
in Section §2.1, the transition in slope of /j, i— > V[fi] along the power curve for 
minima-centered solitons signals a transition from the unstable (positive slope) to 
stable (negative slope) regime. Maxima centered solitons, as discussed, are unstable 
and the transition in slope signals a change in the number of unstable eigenvalues of 
the linearized problem [32, 27]. 

6.3. Numerical methods. The computation of the asymptotic bound states and the "ac- 
tual" bound states are carried out using Matlab and Octave. 

Computation of the Bloch mode, inverse effective mass and . The Bloch mode at 
the band edge, w, is computed using an eigenvalue solver within a single lattice cell (see 
[53] [Appendix] on using Matlab's eigenvalue solver). For convenience we normalize the 
Bloch mode to have unit mass, i.e. J w 2 = 1 . The inverse effective mass tensor A 1 ^ is 
computed by employing Matlab's linear system solver for L" 1 . It is then straightforward to 
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FIG. 6.3. Power— fi plot for Eq. (6.3) with Vq = 10 and K = 2n for bound states centered on a maximum 
(red, dashes) and minimum (blue, solid) of the potential. (Al): Wide view of the semi-infinite gap (semi-log fi- 
axis). (A2) and (A3): Zooming in near the band edge. The asymptotically computed value at the band edge, 
^edge = C* x "Per [Eq. (3.14), black, solid line) and critical power Vcr for the homogeneous (translation invariant) 
equation (dots, black) are shown as well. Geometric shapes correspond to the cases whose bound states are depicted 
in Fig. 6.1. 



compute the inverse effective mass (curvature) 1 , coupling constant 7 e ff, and the band- 
edge power factor using Eqs. (3.7) and (3.15). 

Computation of the bound state at the band edge. The asymptotic bound state is 
comprised of the Bloch mode and the rescaled homogeneous solution. The Bloch mode is 
obtained by periodically extending w from one lattice cell to the domain over which the bound 
state is computed - typically several hundred lattice cells. The rescaled homogeneous ground 
state, F, and its power V C r are computed in ID using the explicit solution, i.e. Eq. (1.13 ) 
with the rescaling in Eq. (6.6). Finally, the asymptotic bound state is obtained by shifting 
F to be centered at point of symmetry of the potential and taking its product with the Bloch 
mode. 

Computation of the "actual" bound states. The bound states of Eq. (1.6) are com- 
puted using Renormalization method [3]. This method is based on fixed-point iterations 
coupled to an algebraic condition, whose role is to constraint the solution to a suitable in- 
tegral identity consistent with the bound-state (otherwise, the iterative solution would di- 
verge). The convergence is monitored by the norm of successive iterations and by 
relative change of the Renormalization constant. For example, for a ID computation with 
|Q| = 0.01 the domain size is a few hundred lattice cells. We use 2 16 grid points to well- 
resolve the oscillations on the scale of the potential period. The computation of the bound 
state is considered to have converged when the difference between successive iterations satis- 
fies |jit n+1 (x) — u™(x)|j oc < 10~ 8 . This typically happens within fewer than 100 iterations 
(a few minutes). 

The Renormalization method needs to be seeded with an initial guess. Deep inside the 
gap the Renormalization method converges when seeded by a Gaussian (or sech) profile. On 
the other hand, near the band edge the method diverges when seeded by a Gaussian or sech, 
which are apparently too far from the basin of attraction of the bound state. We overcome 
this difficulty by seeding the Renormalization method with the asymptotic solution. 

7. Summary and discussion. In this paper we have studied the bifurcation of small 
amplitude (H s (R d ), s > d/2) nonlinear bound states (solitary waves or "solitons") of the 
nonlinear Schrodinger / Gross-Pitaevskii equation with a periodic and symmetric potential. 
Our results provide insight into questions (Q1-Q3) of the introduction. We now briefly sum- 
marize our results, with reference to (Q1-Q3). 
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Concerning (Ql): 

1 . A family of bifurcating solitons (spatially localized standing wave states) can be 
constructed centered at any point of symmetry, xq, of V(x). 

2. Solitons with frequencies near a spectral band edge have a two-scale structure: 
it e (x) w F (e(x — xo)) w(x), where e 2 = — n\ is the distance of the 
frequency to the spectral band edge. 

Concerning (Q2): 

1 . We prove, in general, that the limit of the soliton power, along any family of solitons 
centered at a point of symmetry of V (x), is strictly less than the power of the Townes 
soliton: 

lim V[u(-,n)] = (*Vcr < V cr . (7.1) 

Note: This limit is independent of the centering of the soliton, x . 

2. We prove a high order expansion, which is necessary to capture information about 
the slope of the curve, [i n- V[u(-, near the band edge. Encoded in the slope 
of this curve is information on nonlinear dynamic stability. We conjecture that for 
critical nonlinearities (cr = 2/ d), the curve has positive slope near the band edge 
therefore solitons with frequencies near the band edge are unstable. We have verified 
this analytically for low contrast potentials and numerically for a range of potentials, 
without a smallness constraint on the contrast. 

3. Our analytical results concerning the multiple scale structure of solitons of NLS / 
GP and the curve /i H> "P[it(-, /i)] are corroborated through careful numerical exper- 
iments. 

Concerning (Q3), see Remark 3.3. In particular, see figure 1.1 and the soliton excitation 
threshold conjecture. 

We conclude this section with a discussion the emergent parameter, £*, appearing in 
(7.1). From equations (3.6) and (3.15) we have 



in. 



j B w 2a+2 j V m *y V7cff j voi(B) 



,.,ivf(i£r.fiv/iv > . ( , 2) 



Here m* denotes the determinant of the effective mass tensor, w(x) the B- periodic Bloch 
(band edge) state, 7 e ff the effective nonlinear coupling and vol(B), the volume of the funda- 
mental periodic cell. 

A matter of practical / experimental interest is that the parameters 7 e ff and m* are tunable 

via appropriate design of periodic structure, V(x), therefore making it possible to manipulate 

_ i 

the power curve, V vs. /j,. Figure 7.1 displays m* 2 , and as functions of the potential 
contrast Vq in ID and 2D. All three quantities are bounded between and 1 and decrease 
monotonically with Vq. In particular, this means that Pedge/^cr decreases with Vq. This 
decrease is "faster" in ID than in 2D, at least for Vq < 40. 
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_ 1 

FIG. 7.1. Power factor (blue, solid) and 2 (red, dashes), as functions of the potential contrast Vo in 
(A) ID with cr = 2 [Eg. (6.3) with K = 2n] and (B) 2D with a = 1 [Eq. (6.5) with K x = K y = 2n]. 
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Appendix A. Effective mass tensor. 

In this section we prove equation (3.5), which relates D 2 Ei(Q), the Hessian matrix of 
the band dispersion function E\ to the matrix arising in the multiple scale analysis 

Denote by e l 0(x; k) the Bloch state, associated with E\ (k) : k £ B* —> R; see 
Section §2.2. Here, 0(x; k) is periodic. Thus, 

(- A - 2ik ■ V + |k| 2 + y) 4> = E 1 (k)4> (A.l) 
0(x + q i; k) = 0(x;k), j = l,...,d 

At the bandgap edge one has 

E l (0) = E, , ^(x; 0) =w(x), (A.2) 

where io(x) is the ground state of — A + V, subject to periodic boundary conditions on B. 

We denote . = 9^ . / and f >r = d Xr f. Differentiation of Eq. (A. 1) with respect to ki 
gives 

(-A - 2ihd Xi + |k| 2 + V) (j) kl = (A.3) 

Ei, ki <t> + Ei<f> M + 2iqb, Xi - 2ki<j> , 
At k = 0, Eqs.(A.2) and (A.3) yield 

L^, fci (x;0) = E, M w + 2iw, t . (A.4) 
Removing secular terms from Eq. (A. 6) leads to 



E*m = 0. 
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(A.5) 



It follows from Eqs. (A.4) and (A.5) that 

^ fci (x;0) = 2^L; 1 w ^ . (A.6) 
Differentiating Eq. (A.3) with respect to kj, setting k = and using Eq. (A. 2) we arrive at 

i*0,fc,fc 3 (x;O) = (E, Mkj ~ 25ij) w + [£*, fc >, fej (x;0)+ 

E., kj <f> M teO)] +2i[0 >a:i% (x;O)+0 >% . fci (x;O)] . 

Using Eqs. (A.5) and (A.6) gives 

L*<j) Mk .(x;0) = (E* >k . kj - 2Sij) w-i (d Xi L~ 1 d Xj w + d Xj L~ 1 d Xt w) . 

Removing the secular growth requires that the inner product of the with iu(x) vanish, i.e. 

(E*,kik } ~ 2%) (w,w) - 4:(d Xz L- 1 d Xj w 7 w) 

-4(d Xj L- 1 d Xi w,w) = 0. 

Using the fact that L^ 1 is self-adjoint, the last two terms are equal to each other. Therefore, 
using integration by parts leads to 



1 d 2 E l 



(L^ 1 d x -w.d x w) 



(w, w) 



2 OkiOkj ^ 

This proves the relation (3.5). □ 

Appendix B. Bound on determinant of effective mass tensor. 

Proposition B.l. m^ 1 = dct ( 2 _1 Ei^ kikj (0) ) < 1, with m* = 1 only ifV(x) is 
constant. 

For the proof we use m* > 0; see [34]. Recall that 

1 _ 4(L~ 1 w_j, w.i) 

-E* kiki = Sij — Bij, Bij — -r — ■ — . (B.l) 

2 3 (w,w) 

We claim that B^ is positive definite. To see this, first recall that > with one dimen- 
sional L 2 (T d ) kernel spanned by w. Clearly, w _L M = span {w ; i : i = 1, . . . , d} and 
therefore B is well-defined. 

Let v = («i , . . . , Vd) £ K d be arbitrary. Then 

v • Bv = ( L~\ ■ Vw, v • Vw ) > A^ 1 || v ■ Vw || 2 > C ||vf 

where A2 > denotes the second eigenvalue of acting on L 2 (T d ). 

The matrix E* >kikj is positive definite [34]. Therefore, E* ikikj can be diagonalized by a 
unitary transformation pij such that 

PriE* tki kjPlj = 2 (1 — f3 r ) 5 r i , 
where A-s (i = 1 . . . d) are the eigenvalues of B^. It follows that 
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) = nf =1 (1 - Pi) , 



where fa > 0, i = 1, . . . , d. In order to bound m" 1 from above we will show 
fa € (0,1), i = 1, . . . , d and therefore m" 1 < 1. 

We argue by continuity. Consider the one-parameter family of potentials V(x; 9) = 
9V(x), where 9 G [0, 1] and the associated self-adjoint operator and matrix Bij(9). Since 
L® is self-adjoint and, iu(x; 9), its ground state is simple, there are d continuous functions 
9 fa(6), i = l,...,d, defining the eigenvalues of Bij(9). For 9 = (homogeneous 
medium), E* = 0, Ei(k) = k 2 , and E+fak- — 2%- Therefore B^ = and (3j(0) = 
0, j = l,...,d. In this case (and only in this case!) m^ 1 = dot / = 1. Next consider 
9 = 1, i.e. the original problem. We claim that fail) < 1, i = l,...,d. Otherwise, at 
some value of 9 = 9* > an eigenvalue of B(0*) would attain the value one. This would 
contradict the positive definiteness of Ef k . k .. 

Appendix C. Effective mass for d = 1 and the Floquet-Hill discriminant. In one 

space dimension the endpoints of the spectral bands are obtained by studying the periodic 
and anti -periodic eigenvalue problems [23]. Very briefly, for each E one constructs a 2 x 2 
fundamental solution matrix, M(x; E), and considers the values of E for which M(q\ E) 
has an eigenvalue +1 or —1, corresponding to periodic or antiperiodic eigenvalues. This is 
equivalent to A(E') = ±2, where 

A(E) = trace (M(q; E) = 2cos[k(E)q] . (C.l) 

is the Floquet discriminant. 

The band edge, E = E*, corresponds to k = 0, at which we have A(E*) = 2. Expand- 
ing Eq. (C.l) in Taylor series around k = and E = E* gives 

A(£*) + A' - E.) + 0[(E - E,) 2 ] = 2 ( 1 - + 0(k 4 ) . 



Using A(E^ ) = 2 and solving for the second term on the LHS gives to leading order 



which yields the relation 



p„2 



™- 1 = K(o) = --l q2 



A'(E,) ' 

Since A'(E^ ) < 0, m* > 0. More generally, we have that E"(0) > at the left edge of each 
band and E'-(0) < at the right edge of each band. □ 

Appendix D. Power and slope for small potentials. 

In this section we use a regular perturbation expansion to derive the power and slope 
constants near the band edge, i.e. and assuming a small potential. Such an expansion 
can be made rigorous by an argument based on the implicit function theorem. The derivation 
is comprised of preliminary calculations in any dimension of the Bloch function , an inverse 
linear operator, and the inverse effective mass tensor and coupling constant. To simplify 
notation, subsequent calculations are carried out explicitly in the critical case (d = 1, a = 2). 

REMARK D. 1 . In the derivation below 5 is assumed to be a small constant independently 
of the of t. The calculations are done to order (D(e 2 5 m ) for suitable m. For convenience, the 
e 2 is suppressed from O(-). 
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Let V(x) = 8Vi(x), where \5\ < 1. W.l.o.g. we assume that (Vy) = J B Vi(x) dx = 
Let ws = (x; S), E* — E* (S) be the ground state eigenpair of 



L s w*(x;S) = , L s = L + 6V 1 - E,(S) , L = —A . (D.l) 

with B periodic boundary conditions. 

We expand Wj(x) and E* (<5) in a Taylor expansion in S: 

w s = w + Sw-l + 8 2 w 2 + • ■ ■ , E t (5) = 5E X + 5 2 E 2 + ... , 

where w^ = Wk (x) and we set £^o = since we are interested in bifurcation from the lowest 
band edge. The first three terms in the hierarchy are 



O{6 ) 
OiS 1 ) 

o(s 2 ) 



L Q w Q - , (D.2) 
L oWl = {Ey-V^wo , (D.3) 
L w 2 = E 2 w Q + {Ei - Vi)wi . (D.4) 



Corresponding to the lowest band edge Eq. (D.2) admits a constant solution wo(x) = const 
that spans the kernel of Lq. W.l.o.g. we may choose this constant such that (u>o(x))=l. In 
order to remove secular growth the non-homogeneous terms in Eqs. (D.3) and (D.4) must be 
orthogonal to wo(x). Therefore, their cell-average must vanish. Removing secular terms at 

0(5 1 ) gives 

E 1 = {V 1 ) = 0, w 1 = -Lq 1 V 1 . 
Substituting the above results into Eq. (D.4) and removing secular growth yields 

E 2 = -{VyL^Vy) , W 2 = {L^iVyL^Vy)} , 

where the curly-bracket is a projection symbol defined as 

{/} = /-</>. 

Summarizing the above results gives 

L s = L Q + 5V X + 0(S 2 ) , (D.5) 
w s = l- SL^Vi + 5 2 {L^ l {V l L^ 1 V l ] + 0(S 3 ) , (D.6) 
E*(5) = -^{ViL^Vi) + 0(6 4 ) . (D.l) 

We now approximate the inverse operator LJ . It is expedient to make the following 
definitions. 

DEFINITION D.l. We denote the domains of Lq 1 and LJ 1 as 

Kq = {/(x)|x € B, f is periodic in B, f € Kcr^Lo} , 
Ks = {/(x)|x € B, f is periodic in B, f £ Ker Ls} , 

respectively. We denote the projection operator into Ks as 

D 77. _ 771 (FS,W S ) 

P S F S = F s - rw s . (D.8) 

(w s , w s ) 

Using the above definitions we obtain 
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Lemma D.2. Let F$ = F + SF^ +0{S 2 ). Then 
Lj'PsFg = L^{F } + S [Xq H-Pi} + (WtT^ - Lo^Lo^Fo}}] + 0(6 2 ) 



Proof. The proof of Lemma D.2 follows directly from expanding the projection opera- 
tor (D.8) in powers of 5, using Eq. (D.6), Lq 1 : K -> K , and (L„ V) = 0. □ 

In the derivations of and in each and every case that LJ 1 is applied, one has 
Fo = and ( F±) = 0. Hence, we shall use 

Corollary D.3. It follows from Lemma D.2 that 

F = 0, (Fi) = =>■ L^Fs = SL^Fi + 0(5 2 ) . (D.9) 

In addition, the following approximations of the expressions related to the inverse effective 
mass tensor and effective nonlinear coupling constant are used in the approximation of d*, 

X«(x) = (S* + W, L. ] <), - A i] ) w(x) (D.10) 

= 4 d X] L, l d Xi + x 3 ' ' w x 

\ {w,w) J 

(D.8) , ! v (D.6), (D.9) 2 2 

= Fs (4a Ki L s d x .w§) = -ASd Xi d Xj L Vi + 0{S Z ) , 

„W 1 + 0X11) 



and therefore 



X 2 ( X )^w 2a+1 - lci{ w {D - 6) '^ ,A1} 2aSL^V 1 +0(S 2 ) . (D.12) 



We proceed to calculate the power and its slope. For simplicity we consider the critical 
case (d = l,a = 2). 

Calculation of the power constant . 

Expanding, using Eq. (D.6) with (d = 1, a = 2), gives 



1-25L : 1 V 1 +25 2 {L - 1 (V 1 L Q - 1 V 1 } + 5 2 (L : 1 V 1 ) 2 + 0(6 6 ) 
1 - 6SLq Vi + 65 2 {Lo ^ViLq Vi} + 156 2 (Lo Vi) 2 + 0(<5 3 ) . 



When integrating these functions the contributions of the second and third terms vanish, as 
they are in Kq. Therefore, the first factor in can be approximated by 



Similarly, the inverse effective mass (Gaussian curvature) is 

(d x W,L- l d x w) (D.6) / r-l T/V w r-2 TlM , ^/rfjN 

1-4- r = 1 - 4f (-6d x L Vi)(-5d x L Vi) + 0{6 a ) 



m, 1 = i - 



10, ui) 

2 



1-4^4 d x L^ Vl 
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where in the last step we used the self-adjointness of L . The operator in the above integral 
can be simplified as d x L 3 ^ 2 = Lq 1 . Thus, the second factor in can be approximated 
with 

-L^l-2^/(ioV 1 ) 2 . 
V 71 * jb 

Combining with Eq. (D.13) yields Eq. (3.17). 

Calculation of the slope constant Ci* • 

In one dimension Eq. (3.16) reduces to 



Ci* = 



/ / \U 1 (x,y)\ 2 dxdy + -f w 2 dx f S(y) d u F(y; 0)| fi= _! dy . (D.14) 
JrJb Jb Jm. 



IRJB 

(1.13) 



As we shall see, U\ = 0(S) and S(y) = 0(5 2 ). Therefore, both integral terms are 0(8 2 ). 
It follows from Eqs. (3.4), (D.10), and (D.ll) that to leading order in 8, F(y) is the Townes 
mode, i.e. 

F(y;Sl,S) = R(y;n) + 0(5 2 ) . (DAS) 

Furthermore, we shall use Eq. (1.13) with a = 2 to explicitly evaluate the y-integrals. 
The first integral term (D. 14) depends on 

Ui(x,y) {3 = 1] 2Lj 1 d x wd y F ( °= 6) ^^{Sd^V^Ry + 0(S 2 ) 

(I ^ 9) -2S8 x Lq 2 ViR y + 0(5 2 ) . (D.16) 

Therefore, 

/ / If/i^y)! 2 dxdy = \8 2 l {d^V^ 2 dx f R 2 y dy + 0(8 3 ) 

JRJB JB JR 

V3 tt<5 2 / {d x L^ 2 Vif dx + 0(5 3 ) . 
Jb 

The operator in the above integral can be simplified as d x L^ 2 = (— d xxx ) _1 . Hence, 

/ / ^{x^ 2 dxdy = V3tt5 2 J [(-d^r 1 ^] 2 dx + 0(6 3 ) . (D.17) 

JRJB JB 

For the second integral term in (D.14) we need to calculate S(y). We use 

U (x, y) ( = 6) w(x)F(y) = R(y) + 0(S 2 ). (D.18) 

The first two terms in Eq. (3.13) are negligible. This follows from 
Lemma D.4. (U 2p , w) = . 

Proof. Equation (4. 16) shows that the .x-dependence of U-2 P is of the form L~ 1 Xk(x) for 
suitable Xk{x), k = 1, 2. As LJ 1 is into the orthogonal space to w(x), the Lemma follows. 
□ 

That the first term in Eq. (3.13) is zero follows immediately from Lemma D.4. The sec- 
ond term in Eq. (3.13) has an additional C/q' 7 . However in light of Eq. (D.18) the x— dependence 
of Uq 17 is constant to leading order. Therefore, the second term in S(y) is 0(6 3 ). It remains 
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to calculate the two last terms in S(y). Note that the coefficient preceding the square brackets 
in Eq. (3.13) cancels with the iw-integral in Eq. (D.14). This leaves (assuming d = 1, a = 2) 

S(y) = W(w,U$U 2 ) - 2^d x w,d y U^ + 0(S 2 ) , (D.19) 



where 



9 X Lq 1 XlOyyyF + 8 X Lq 1 X 2 Oy F 

+Lo L d x w(d vy ~ l)d y F + 2w 4 Lo 1 d x wd y F* 

and in one-dimension Xi = X^ . We denote by Sk the various terms in S(y) when C/3 is 
explicitly inserted into it, and by If* their corresponding contributions to Eq. (D.14). 
Using Eqs. (D.6), (D.16) and (D.18), the first term in Eq. (D.19) is 

St = 10(w, UgU?) = lQ(Uf) + 0(5 3 ) . 

Therefore, its contribution to the slope is 10 times the first term in Eq. (D.14), i.e., 

7j (D = 17) 10V3 n8 2 -f [(-Qbxx) -1 ^] 2 dx + 0(5 3 ) . (D.20) 
Jb 

The first term arising from substituting t/3 into S(y) is 
S 2 = -2(d x w,d y 2LQ 1 d x L^ 1 X 1 dy yyy F) 

(D.6)^D.10) 2 (_ §L -i dx y u2L -i dxL -i (-ASd^L^V^dyyyyF + 0(S 3 ) 



( D J 5 ) .axlf fx r-5/2 



2 



lBSy [d XX L °^Vlj dxdyyyyR+0(S 3 ) , 

where in the last step we used the skew self-adjointness of d x and the self-adjointness and 
positivity of Lq . Simplifying the operator in the above integral we obtain 7 

S 2 = -l^l [(-d^yWt] 2 dxdyyyyR+0(S 3 ) . 

Jb 

Substituting S 2 into Eq. (D.14) gives 

h = -m 2 i [(-d^r 1 ^] 2 dx [ dyyyyR d n R{y;£l)\n=~i dy + 0{5 3 ) . 
Jb Jr 

The following explicit integral can be obtained from Eq. (1.13) 



dyyy y R d n R(y,Q)\ n= _ x dy 



16 



Using this gives 



h = nVsnS 2 ! [(-a^rVi] 2 dx + o(6 3 ) . (D.21) 

Jb 



7 2 5 /2 

'Note that the Fourier representations of d x L Q andd xx L are the same only in d = 1. 
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Similar calculations can be carried out for I 3 , 14 and I 5 using the explicit integral 

/ d yy R 5 d Q R(y;Q)\ a= _ 1 dy=^^. 

Thus, to 0(S 3 ), we get 

h = 13x/3^ 2 / [(-d^)^] 2 dx 

JB 

and ^5 = —I4 = jl 3 . Summing the contributions from I\ . . . I§ gives 

Ci* = 34V3^ 2 / [{-d xxx )- x V^ dx + 0{5 3 ) > . 

JB 

This concludes the proof of Corollary 3.3. 
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